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ABSTRACT 


Progress  toward  the  development  of  an  electromagnetic  pulse 
sounding  probe  for  the  detection  and  delineation  of  geological 
anomalies,  primarily  in  a  hard  rock  medium,  is  reported.  This 
report  covers  research  performed  on  Contract  H0210042  during  the 
period  24  February  1971  to  23  February  1972.  Analytical  studies 
on  design  data  for  the  probe  and  on  the  scattering  by  planar  and 
spherical  contrasts  are  summarized.  Successful  measurements  made 
using  a  "first  generation"  version  of  the  probe  on  conducting 
and  dielectric  cylindrical  targets  in  a  lossy  overburden  and  on 
specific  anomalies  in  a  limestone  medium  are  reported. 
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I.  TECHNICAL  REPORT  SUMMARY 

diagnosis*and  5r0be  for  the  Action, 

exploration]  geophysical  "prospect  ingSandhas  appl^f‘ ions  ^^geophysical 

and  interrogating  pulse  wild  be  aSSsteTtS  m  fh£r°be  St[UCtUre 
the  problem  The  anaic  nf  +hQ  ajustea  to  fit  the  parameters  of 

are  to  design  and  build  a  nrnhpP!:t9ra+1  dun!??  this  contract  period 

electromagnetic  energy  into  a  arnn^C*Ure  “Jich  effectively  couples 

<stics9„fsejected^^ra^nri;i^r-j^'b--r. 

two  partS?ffFirst!ranalytical0studipsPwe10d  ta5  been  separated  f"‘° 
data  for  the  probe  and  to  obtain  thpSQp»ttpC™dUCJed  t0  secure  design 
planar  and  spherical  contrast  L  -a  s??ttdn"9  characteristics  of 
anomalies.  Jt  ti  a  !  '  t  'deal, zed  models  of  geological 

and  interfacing  to  initiate  an  experiment!!  ePu'Pment,  instrumentation 
assembled  and  tested  The  spcnnr/mJ?6^3  measurement  program  were 
ment  program,  was  viewed  as  thp  mnEf  °f  ?^r  effort,  the  measure- 
tract.  Our  reasoning  was  thaJ  t^SJ  ***$}'•?  portion  of  the  c°"- 
technique  must  be  demonstrated  first^in^  Veal°w,  PUl-6  SOUndl’ng 

using  admittedly  less  than  nntim,™  a  ?  1  world  environment, 

methods.  Once  this  was  dSnpP  an,7wf-0bf  d£Tgns  a,?d  Pressing 
optimize  the  performance  of  thp  nuiV03  aJ]d  exPen’mental  studies  to 
on  later  contracts  with  the  securp  knnwipHdinLPr°be  could  be  Pursued 
tainable.  For  this  reason  nn r*  t  !^wl?dge  that  the  Seals  are  at- 
probe  could  be  constructed,’  the  con  Jacf  Versi'?n  °f  the 

concentrated  on  the  measurement  effort  Snseo^tlT"^"61  Wer+ 
programs  developed  to  nrnvirip  nmh0  !,*•  consequently,  the  computer 

not  been  fully  exploited  as  "et"  9"  and  SCatten'ng  data  have 

a  source  which  producesCaPperiodicntrainSof ^ tern  currently  in  use  employs 
pulses  are  coupled  into  the  SrounS  nr  In  f  v^eo-type  pulses.  The 
arm  of  a  broadband  probe  consistina  nf  a  k  ?edl,ini  by  the  transmit 

balun-fed  dipoles  stretched  over  the  Surface *"0?  thptIlpH°nallyT?riented 
probe  can  be  operated  either  in  *  5  uthe  medlum-  The 

IS  used  or  in  aT ort^^e on^n^  °",y  one  dip0'e 
and  the  other  as  a  receiver  In  the  nrthLn  6  serves  as  a  source 
not  respond  to  the  air-medium  infprfanpth°;?0na1  mode’  the  Pnobe  does 

surface  obstructions  (a  vehicle 
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not  detected).  On  a  rock  medium,  some  sensitivity  to  large  surface 
obstructions  is  noted,  but  a  change  in  the  shape  of  the  dipo  e  arms 
should  obviate  this  problem.  In  operation,  a  pulse  is  launched  on 
one  dipole  at  the  feed  point  which  penetrates  into  the  medium.  As 
the  wave  propagates  along  the  dipole,  its  fringe  fields  penetrate 
deeper  into  the  medium.  Reflected  signals  from  the  medium  are 
observed  on  a  sampling  scope  either  on  the  same  dipole  (direct 
mode)  or  on  the  second  dipole  (orthogonal  mode)  acting  as  a 
receiver.  In  the  orthogonal  mode,  the  probe  is  "blind"  to  ho¬ 
mogeneous  targets,  e.g.,  horizontal  stratifications,  whereas  in 
the  direct  mode  both  homogeneous  and  inhomogeneous  targets  are 
seen. 


All  measurements  reported  for  this  period  were  made  with  a 
5  volt  peak  pulse  generator  (3  ns  base,  2  MHz  repetition  rate). 

The  fact  that  targets  could  be  detected  in  both  the  overburden 
(depths  to  5  feet)  and  hard  rock  (depths  to  10  feet)  using  such 
low  power  is  extremely  encouraginq.  New  pulse  generators*  and  better 
baluns  now  available  at  the  ElectroScience  Laboratory  will  give  us 
both  higher  power  and  lower  frequencies  in  the  spectrum  (lower 
repetition  rates);  we  anticipate  successful  target  detection  at 
significantly  greater  depths  with  this  equipment. 

fTC:  mquybtisgp  the  -f9chnic«i  Project  Officer,  the 

experiments  discussed  herein  are  directed  toward  the  detection  and 
identification  of  buried  cylinders.  In  addition  to  numerous  probe  tests 
and  control  or  no-target  measurements,  the  following  targets  have  been 
measured  using  the  probe  in  the  orthogonal  mode.  1,  3,  5  and  10  foot 
long  hollow  metal  cylinders  with  diameters  of  4  inches  buried  at  a  depth 
of  5  feet  in  the  overburden.  2,  3,  and  4  inch  diameter  extended  plastic 
pipe  (hollow)  buried  at  a  depth  of  3  feet  in  the  overburden.  A  "dummy" 
trench,  i.e.,  a  trench  identical  to  those  dug  to  install  the  cylinder 
targets  was  dug  (5  foot  depth)  and  then  refilled.  In  all  cases,  the 
refilling  was  compacted.  At  a  local  limestone  quarry,  measurements 
were  made  of  a  hollow  metal  pipe,  1  foot  in  diameter  whose  top  was 
buried  2.0  feet  in  a  limestone  medium,  inhomogeneous  stratifications 
in  a  limestone  medium  and  a  large  tunnel  (road)  beneath  some  30  feet 
of  limestone.  Several  control  measurements  were  also  made  for  the 
limestone  medium.  On  May  31,  1972,  after  this  research  and  this  report 
had  been  completed  the  authors  were  informed  that  the  interest  of  the 
Technical  Project  Officer  had  shifted  so  that  extended  planar  targets 
were  of  major  interest.  An  immediate  conclusion  can  be  made,  on  the 
basis  of  scattering  theory,  that  the  scattered  pulses  for  extended 
planar  targets  should  be  substantially  larger  than  those  from  the 
cylinders  treated  in  this  report. 

For  both  the  limestone  and  overburden  measurements,  the  target 
was  detected  in  each  case  with  the  exception  of  the  tunnel.  We 
would  not  realistically  expect  to  penetrate  some  60  feet  (2-way  path) 
of  limestone  with  only  a  5  volt  peak  pulse  with  frequency  content 
above  1  MHz.  The  presence  of  the  tunnel  should  be  easily  evident 
with  measurements  using  high  prwer  pulse  generators.  Presence  of 


*A  Hewlitt  Packard  model  2/4A  pulse  generator  was  purchased  with 
contract  funds. 
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th“  target  was  evident  both  from  a  simple  comparison  of  the  tarqet 
and  control  time  domain  waveforms  and  from  examination  of  the  spectral 
components  of  these  waveforms.  In  addition,  the  reflected  signals  from 
the  pipe  and  inhomogeneous  stratification  targets  in  limestone  and  the 
pla  tic  pipe,  cylinders  and  dummy  trench  were  separable  in  the  sense  that 
the  scattered  waveforms  and  their  associated  spectra  contained  adequate 
information  to  identify  the  individual  targets  in  the  set  of  targets 
However,  with  the  simple  processing  used  it  was  not  conclusively 
evident  that  one  could  separate  the  different  lengths  of  cylinder 
This  is  expected  since  a  theoretical  analysis  shows  that  the  current 
distribution  on  a  buried  wire  antenna  is  of  a  non-resonant  form 
Thus  the  only  identifiable  characteristic  of  the  linear  scatterers 
would  be  the  relative  magnitude  of  the  scattered  field  and  this  is 
a  very  poor  characteristic  for  the  purposes  of  identification. 

Since  the  cylinder  diameter  is  quite  small,  less  than  4  inches, 
separation  of  various  cylinder  diameters  is  not  important  for  the 
purposes  of  this  study.  Measurements  of  a  larger  cylinder  at  the 
limestone  quarry  do  indeed  indicate  that  separation  is  possible  in 
terms  of  cylinder  diameter. 

HTl!e  ana]/tic^  tools  for  the  design  of  an  electromagnetic  pulse 
sounding  probe  and  for  the  analysis  of  planar  and  spherical  sub¬ 
surface  targets  have  been  developed.  Extensive  computations  usinq 
computer  programs  developed  to  exploit  these  analytical  tools  have 
been  deferred,  however,  in  order  to  concentrate  on  an  initial  series 
of  experimental  measurements.  A  "first  generation"  version  of  the 
pulse  sounding  probe  has  been  designed  and  built.  A  mobile  riq 
which  permits  measurements  to  be  made  at  remote  sites  (removed  from 

fMrhichoIiteriia"d  w^ho^.avai1ab^e  electrical  power)  has  been  re¬ 
furbished.  Using  the  first  generation"  probe  and  a  low  power 

pulse  generator,  successful  detection  of  conducting  and  dielectric 
cylinders  at  depths  to  5  feet  in  the  overburden  and  of  a  metallic 
cylinder  and  inhomogeneous  stratifications  to  depths  of  10  feet  in 
a  limestone  medium  has  been  demonstrated.  The  necessary  equipment 
o  operate  the  probe  at  significantly  higher  power  with  pulse  gen¬ 
erators  having  much  lower  spectral  content  has  been  acquired  and  is 
ready  for  testing.  The  feasibility  of  an  electromagnetic  pulse 
sounuing  technique  has  been  established  and  we  are  now  ready  to 
refine  our  methods  and  equipment  for  specific  applications  to  the 
hard  rock  tunneling  problem. 


II.  PURPOSE 

j •  A  ^nicIu®  electromagnetic  pulse  sounding  probe  for  the  detection, 
J  n°sl  and  identification  of  geological  and  man-made  anomalies 

thin  the  earth  is  being  developed.  The  immediate  application  of 
the  probe  is  as  a  hazard  detection  device  in  advance  of  hard  rock 
rapid  tunneling  operations,  but  numerous  other  applications  in 
geophysical  exploration,  geophysical  prospecting  and  in  mininq  or 
dnlling  operations  can  be  envisioned.  During  this  contract  period 

neh°bj?ctl!es  of  the  program  are  to;  1)  design  and  demonstrate  a 
probe  structure  capable  of  effectively  coupling  electromagnetic 
pulse  energy  into  a  ground  or  rock  medium  and  at  the  same' time  be 
insensitive  to  obstructions  located  on  the  medium  surface. 


2)  interrogate  with  this  probe  certain  simple  targets  in  ground  and 
rock  media  and  3)  establish  preliminary  bounds  on  the  workable  depths 
of  electromagnetic  pulse  sounding  via  theoretical  and  experimental 
methods. 


III.  INTRODUCTION 

The  ultimate  goal  of  the  research  on  this  program  is  to  in¬ 
corporate  state-of-the-art  advances  developed  at  this  laboratory, 
both  in  antenna  systems  and  processing  techniques,  into  a  geological 
tool  capable  of  delineating  the  size  and  shape  of  geological  anomalies 
from  the  shape  and  spectral  content  of  an  echo  pulse. 

The  detection  of  underground  anomalies  that  represent  hazards 
using  a  periodic  video  pulse  radar  system  does  not  appear  to 
present  substantial  difficulties.  The  problems  that  are  of  most 
importance  are  those  associated  with  isolating  the  scattered  pulse 
associated  with  buried  target  and  then  its  identification.  Pulses 
can  be  scattered  from  many  different  discontinuities  including 
scattering  from  the  balun,  the  input  terminals  of  the  antenna,  the 
ground  surface,  targets  above  ground,  etc.  Much  of  our  effort  has  been 
focussed  on  reducing  these  undesired  scattered  pulses  so  that  the 
scattered  pulses  from  the  underground  target  (lizard)  can  be  isolated. 
This  effort  has  proven  successful.  Once  the  scattered  pulse  from 
the  underground  is  properly  separated  from  the  clutter,  there  still 
remains  the  problem  of  identification.  While  only  a  limited  effort 
has  been  expended  for  this  purpose  under  the  current  contract,  there 
has  been  developed  at  the  ElectroScience  Laboratory  a  variety  of 
techniques  for  target  identification.  Since  these  techniques  are 
required  in  a  successful  system,  they  have  strongly  influenced  the 
choice  of  parameters  of  the  system  to  be  used  for  the  electromagnetic 
hazard  detection  system. 

Studies  at  this  laboratory  have  conclusively  demonstrated  that 
if  the  gross  features,  i.e.,  the  overall  size  and  shape,  of  a  target 
are  to  be  deduced  via  electromagnetic  interrogation  then  discrete 
frequencies  over  an  8  or  10:1  bandwidth  with  a  fundamental  such  that 
the  linear  extent  of  the  anomaly  is  approximately  1/10  of  a  wavelength 
in  the  ambient  medium  are  dictated J  With  a  pulse  sounding  approach* 
to  detection  of  subsurface  targets,  this  fixes  the  characteristics 
(power,  rise  time,  repetetion  rate,  etc.)  of  a  pulse  generator  suitable 
for  interrogation  of  a  given  range  of  anomaly  sizes  at  a  given  depth. 

It  does  not  necessarily  follow  that  all  of  the  higher  spectral  con¬ 
tent  are  needed  to  obtain  a  viable  signature,  but  our  experience 
indicates  that  the  lowest  frequencies  are  essential. 2  Of  course,  the 
more  high  frequency  content  available,  the  more  detail  of  the  target 
is  available,  and  the  sharper  the  reflected  pulses  yielding  an  im¬ 
proved  range  gating.  If  follows  that  more  detail  and  an  improved 


*An  alternative  system  for  radar  applications  uses  discrete 
harmonically  related  cw  measurements  over  the  10:1  bandwidth. 3 


4 


range  gating  a^e  to  be  expected  when  comparing  measurements  in  a  hard 
rock  medium  to  those  in  the  overburden  since  the  high  frequency  energy 
is  severely  attenuated  in  the  overburden.  The  use  of  low  freguencies 
to  interrogate  the  target  is  ideally  suited  to  the  problem  of  sub¬ 
surface  targets  in  either  a  ground  or  rock  medium  since  it  is  precisely 
these  frequencies  which  can  penetrate  to  the  target. 

for  the  case  of  targets  in  free  space,  i.e. ,  radar-type  targets, 
studies  at  this  laboratory  have  established  three  distinct  techniques 
for  processing  low  frequency  scattering  data  to  def“^t  and  identify 
the  target.  One  of  these  uses  harmonic  spectral  amplitude  data  in 
an  n-dimensional  space  approach  to  target  classification. 4  A  second 
approach  uses  harmonic  spectral  amplitude  and  phase  data  to  derive 
a  small  set  of  excitation-invariant  target  descriptors.5  In  the 
third  approach,  relationships  between  the  size  and  shape  of  the 
target  and  its  low  frequency  scattering  characteristics  are  exploited 
to  obtain  3-dimensional  isometric  images  of  the  target  from  harmonic 
spectral  amplitude  and  phase  data.3  Ultimately  it  is  intended  to 
apply  these  and  similar  sophisticated  processing  schemes  to  the 
analysis  of  electromagnetic  pulse  sounding  data.*  At  this  stage, 
however,  we  have  confined  ourselves  to  relatively  simple  averaging 
and  differencing  techniques  with  subsequent  comparisons  of  target 
and  control  measurements  in  both  the  time  and  frequency  domains. 

Once  our  probe  designs  have  been  optimized  for  a  given  medium  and 
techniques  for  handling  high  power  pulses  with  very  broad  spectral 
content  established,  we  can  properly  begin  to  consider  more  ad¬ 
vanced  processing. 

In  Section  IV  of  this  report,  the  status  of  our  analytical 
studies  and  the  availability  of  computational  programs  are  sum¬ 
marized.  Details  of  the  analytical  studies  have  been  given  in  a 
previous  report,  which  is  included  as  an  Appendix  of  this  report. 

Section  V  of  this  report  describes  the  electromagnetic  pulse 
sounding  probe,  reviews  the  probe  design,  tests  and  operation  in 
both  the  direct  and  orthogonal  modes;  and  describes  the  processing 
used  on  the  pulse  sounding  data. 

In  Sections  VI  and  VII  respectively,  experimental  data  on  man¬ 
made  targets  in  a  lossy  overburden  and  man-made  and  natural  targets 
in  a  limestone  medium  are  presented. 

Section  VIII  of  the  report  sunmarizes  our  conclusions  regarding 
the  feasibility  and  applicability  of  an  electromagnetic  pulse 
sounding  system  for  interrogating  subsurface  targets  and  makes 
specific  recommendations  for  future  research  in  this  area. 


*The  n-dimensional  space  approach  has  previously  been  used  with 
shallow  subsurface  target  measurements. 6 
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IV.  THEORETICAL  STUDIES 


Two  types  of  theoretical  studies  have  been  conducted  on  the 
program.  The  first  was  a  comprehensive  analysis  and  computer 
program  for  all  types  of  wire  antennas  and  arrays  of  wire  antennas 
in  a  homogeneous  dissipative  medium.  The  second  type  was  plane 
wave  scattering  computations  for  planar  and  spherical  conductivity 
and  permittivity  contrasts.  In  addition,  closed  form  results  for 
particular  orientations  for  the  transient  fields  of  an  infinite 
line  source  above  a  dielectric  half-space  were  obtained.  Approxi¬ 
mate  results  estimating  the  effects  of  the  overburden  on  the  scat¬ 
tering  from  various  targets  were  obtained.  Illustrative  results 
from  the  above  studies  were  given  in  the  Semi-Annual  Technical 
Report  which  is  included  as  an  Appendix  of  this  report. 

..  .  Ljsted.  bel°”  ar?  brief  descriptions  of  the  computer  programs 
which  have  been  developed  and  are  now  available  to  provide  probe 
design  data  and  scattering  computations. 

1.  Wire  Antennas  in  a  Homogeneous 
Dissipative  Medium 

HinnW6  5r^raT  hand1es  a11  types  of  wire  antennas  including  linear 
dipoles,  V  dipoles,  rectangular  loops,  circular  loops  and  arrays, 
he  wire  antenna  can  have  a  finite  conductivity,  and  certain  portions 
of  the  antenna  can  be  enclosed  by  a  thin  dielectric  sleeve.  The  out¬ 
put  from  the  program  provides  a  complete  analysis  of  the  system  in¬ 
cluding  self  impedance,  mutual  impedance,  current  distributions  near- 
??!J®  ?S  and  far"fie1d  Patterns.  Solutions  are  via  a  piecewise- 
sinusoidal  expansion  for  the  unknown  current  distribution  and 
Galerkin  s  method  for  reducing  the  integral  equation  formulation  to 
a  system  of  simultaneous  linear  equations.  The  present  analysis  and 
programs  require  the  medium  to  be  of  infinite  extent.  For  prow r 
probe  design  data,  these  results  need  to  be  extended  to  incfude  an  air- 
dissipative  medium  interface.  The  first  attempt  to  do  this  via  an 

b^ausrthpr^nf°rmIi1atl0n/lth-1ar9e  interelement  spacing  was  discarded 
antfnf6  reP|jired  computer  time  became  prohibitive.  For  finite 

nrnhlprK  °f  arb^ary  sh?Pe»  the  half-space  introduces  rather  difficult 

fundin9  levei  and  the  „ 


-•  Homogeneous  Dissipative  Spherical  Targpt. 
in  a  Homogeneous  bissir~'*'~~~~  u  J'’~ - — 


i pative  Medium 

f°9ram  calculates  the  scattered  fields  (plane  wave  in- 
cidence)  from  an  arbitrary  homogeneous  spherical  target  with  con- 

™£Pc™oterS  i-er?ed  in  an  -bl2t  JSfl'Sth 

parameters  c2,  o2,  uQ.  Solution  is  via  a  Mie  series  formulation. 

3*  Homogeneous  Layered  Dissipative  Spherical  Target 
__in  a  Homogeneous  Dissipative  Medium - 

?r°9ram  calculates  the  scattered  fields  (plane  wave  in¬ 
cidence)  from  an  arbitrary  homogeneous  spherical  target  (q,o1)lJo) 

G 


witF  a  homogeneous  spherical  layer  (c9,o9,u)  irrmersed  in  an 
horn  leneous  ambient  medium  Uo.Qo.uJ?  §ol8tion  is  via  a  Mi< 
series  formulation.  J  J  0 


4 .  Plana  r_  Co n trasts 

The  program  calculates  the  impulse  response  wavefonn  (inverse 
transform  of  the  frequency-dependent  Fresnel  reflection  coefficient) 
for  normal  incidence  at  the  interface  of  two  homogeneous  media  ,i 
<l,uo  and  r2>°2»np*  The  singularities  and  discontinuities  in  the* 
waveform  can  be  obtained  in  closed  form.  The  program  yields,  via  a 
Fourier  synthesis  procedure,  the  decay  of  the  waveform  from  its 
initial  value.  Frequency-dependent  constitutive  parameters  can  be 

5*  Interface  and  Overburden  Effects 

The  program  calculates  the  back  scattered  field  from  a  sub¬ 
surface  target  assuming  off-normal  plane  wave  incidence  on  the  air- 
ground  interface  and  frequency-independent  constitutive  parameters 
for  the  ground  medium.  It  is  also  assumed  that  the  target  scatters 

niSjtn..t1Hf1JhSCalin9  °f  •itSufree"sPace  scattering  characteristics. 

°V!?e  P^ram  is  the  ramp  response  waveform  (impulse  response 
waveform  twice  integrated  with  respect  to  time)  of  the  target  and  this 

waveform  is  compared  with  the  free  space  ramp  response  waveform 
ot  tne  target, 

6.  Convolution 

The  program  calculates  the  convolution  of  two  time  domain 

S*""™!:,.  wavef0™  completely  arbitrary  and  is  approximated 
by  straight  line  segments  At  present,  the  second  wavefonn  must 
be  twice  integrable  analytically  in  closed  form,  e.g.,  a  pulse,  pulsed 

innrHnn  +The  ™te"ded  use  of  -his  Program  is  primarily  in  con¬ 
junction  with  targets  whose  impulse  or  baseband  pulse  response  is 

known  only  as  a  graph  or  plot,  i.e.,  not  analytically.  Using  the 
quickly'obtained961*5  reSP°nSe  to  other  interrogating  signals  can  be 

^nother  Program  used  to  control  our  pulse  sounding  measurements 
and  to  process  and  manipulate  the  measured  data  is  described  in  a 

withVn^ILfl  f£nn  reporJ*  The  Programs  listed  above  are  used 
a  cmaii0^  J  ft  J6??  comPuter*  The  control  program  is  written  for 
a  small  instrumentation  computer  (IBM  Minimal  Informer)  which  is 

f°r  .on~]ine  corT1Puter  controlled  measurements.  The 
instrumentation  and  software  to  shift  the  on-line  measurements  to  the 
more  versatile  modified  Oatacraft  System  are  being  prepared  at  this 

a.  Probe  Design 

, ,  j*1®  W1re  antenna  analysis  and  programs  provide  a  straightforward 
method  for  optimizing  the  probe  design  for  a  given  target  range 
frequency  range  and  medium  conditions.  These  are  based?  at  the’ 
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fi cations  of  the  program  to  account  for  the  ai '''-medium  interface  are 
not  yet  complete.  Experimentally  we  have  demonstrated  that  the  mis¬ 
match  due  to  the  interface  can  be  controlled  by  adjusting  the  shape 
of  the  dipole  arrs  and  the  depth  of  grounding  stakes.  Once  we  have 
analytically  accounted  for  the  interface,  the  dipole  arm  shape  and 
ground  rod  depth  can  be  optimized  for  a  w’de  variety  of  media. 

The  analyses  and  computer  programs  listed  above  assume  that  the 
conductivity  and  dielactric  constant  are  constant  with  respect  to 
frequency.  It  is  recognized  that  these  parameters  are  frequency 
dependent.  The  introduction  of  complex  frequency  dependent  con¬ 
stitutive  parameters  is  a  simple  matter  of  changing  values  for  frequency 
domain  computations.  Periodic  time  domain  waveforms  can  be  constructed 
from  these  computations  using  the  available  Fast  Fourier  Transform 
(FFT)  computer  program,  if  needed.  This  cannot  be  readily  achieved  for 
the  direct  time  domain  computations  of  the  waveforms.  In  most  in¬ 
stances,  this  is  not  an  essential  flaw,  however,  since  the  major 
variation  is  of  the  form 


where  c* ,  the  complex  dielectric  constant  is  frequency  dependent. 

In  a  previous  report,7  calculations  using  the  wire  antenna 
programs  demonstrated  that  for  shallow  depths  (<  4  feet)  an  un¬ 
insulated  dipole  produced  field  strengths  an  order  of  magnitude 
greater  than  those  of  an  insulated  dipole.  However,  at  much  greater 
depths  only  very  slight  differences  in  field  strength  are  noted. 

In  general,  we  are  interested  in  targets  at  some  distance  and 
conclude  that  for  our  present  applications,  either  insulated  or 
uninsulated  probes  can  be  used. 

In  Fig.  1,  the  near-zone  fields  in  the  plane  of  a  U-shaped 
dipole  in  a  homogeneous  ambient  medium  cr  =  9,  o  =  0.003  mho/meter 
are  shown  at  ranges  of  15  and  30  feet.  The  dipole  has  15  foot  in¬ 
sulated  arms  with  2.5  foot  uninsulated  ends  perpendicular  to  the 
main  arms.  The  results  in  Fig.  1  were  obtained  from  the  wire 
antenna  program  at  a  frequency  of  10  kHz.  It  was  found  that  the 
results  remained  essentially  unchanged  up  to  a  frequency  of  100  kHz. 
Note  that  the  field  component  Ex  which  is  parallel  to  the  dipole 
arms  has  a  sharp  null  which  moves  with  range.  This  would  be  of  no 
consequence  for  extended  targets,  but  for  finite  objects  the  dipole 
would  have  a  "blind  spot"  in  the  direct  reflection  mode  for  the 
x-component  of  the  field.  In  Fig.  2,  the  same  results  are  shown 
where  the  uninsulated  portion  of  the  dipole  arms  have  been  in¬ 
creased  to  5  feet.  Note  that  at  the  15  foot  range  lengthening  the 
uninsulated  portion  of  the  dipole  arms  has  simply  raised  the  level 
of  the  x  and  z  components.  The  same  is  true  at  the  30  foot  range. 

In  each  case  a  3  to  4%  increase  in  field  strength  is  noted. 

The  wire  antenna  program  permits  one  to  calculate,  for  a 
homogeneous  dissipative  medium,  the  coupling  between  two  arbitrary 
wire  antennas  separated  by  a  given  distance  and  in  any  orientation. 
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Fig.  1.  Calculated  near  zone  fields  of  a  wire  antenna 
in  a  homogeneous  conducting  medium. 

(b)  30  foot  range. 
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Fig.  2.  Calculated  near  zone  fields  of  a  wire  antenna 
in  a  homogeneous  conducting  medium. 

(a)  15  foot  range. 


Fig.  2.  Calculated  near  zone  fields  of  a  wire  artenna 
in  a  homogeneous  conducting  medium. 

(b)  30  foot  range. 


The  received  voltage  at  a  second  dipole  when  the  first  dipole  is 
excited  by  a  1  volt  input  is  shown  in  Fins.  3,  4  and  5  over  a  fre¬ 
quency  range  from  10  KHz  to  10  MHz.  The  dipole  configurations  are 
shown  in  the  figures,  and  results  are  shown  for  various  assumed 
permittivities  and  conductivities  for  the  medium.  In  Fig  3  the 
re  at ive  permittivity  is  15.0  and  results  for  conductivities ’of 
0.1,  0.05  and  0.01  mhos/meter  are  given.  The  dipoles  have  10  foot 
arms  and  are  separated  by  approximately  5.4  feet.  In  Fig.  4,  the 
dipole  arms  are  reduced  to  5  feet,  the  relative  permittivity  is  15.0, 
and  the  same  conductivities  are  used.  Also  shown  is  the  result  for 
a  restive  permittivity  of  9.0  and  conductivity  of  0.01  mhos/meter 
llTustrating  that  below  1  MHz  the  results  are  essentially  independent 
o  the  dielectric  constant.  In  Fig.  5,  a  portion  of  one  antenna  is 
insulated  and  the  arms  (10  feet)  are  in  a  U  configuration.  The 

nenn^iVeJPnr^lotlvity  is  9,0  and  resu1ts  for  conductivities  of  0.01, 

0  006  and  0.003  mhos/meter  are  shown.  Note  from  all  three  figures 
that  over  a  lower  portion  of  the  spectrum,  the  received  voltage 
actually  increases  with  frequency  to  some  peak  value  and  then  de¬ 
creases.  Also,  this  effect  becomes  more  pronounced  as  the  con¬ 
ductivity  decreases. 

Oalcinations  such  as  those  shown  in  Figs.  3,  4  and  5  permit  one 
to  predict  the  attenuation  as  a  function  of  frequency  for  a  medium 
with  a  given  dielectric  constant  and  conductivity  by  comparing 
calculations  at  2  different  ranges.  In  Fig.  6,  the  predicted' attenuation 
in  dB  per  foot  is  shown  over  the  frequency  range  10  KHz  to  10  MHz  for 
a  medium  with  a  relative  permittivity  of  15  and  conductivities  of  0.1, 

?An5Lund*?'°  mhos/fneter-  Note  that  over  the  frequency  range  10  kHz  to 
400  kHz  the  curves  are  essentially  flat,  i.e.,  the  attenuation  is 
constant.* 

The  attenuation  of  the  local  ground  has  been  measured  in  situ 
using  pulse  propagation  data  between  a  dipole  buried  at  a  depth  of 
5  feet  in  the  overburden  and  3  dipoles  on  the  surface.  An  averaqe 

romrinekHza  n  ??ta  are  sh0Wn  in  Fi^  7  for  frequencies 

IP  kJz  t0  0  MHz*  Actually,  two  pulse  generators  were  used  to 

span  the  frequency  range  and  the  two  points  shown  for  1  MHz  illustrate 

dataSineFin0f7th-tha^  from  *he  different  generators.  Comparing  the 
data  in  Fig.  7  with  the  predicted  attenuation  curves  in  Fig.  6  it 

is  clear  that  a  qualitative  agreement  between  predicted  and  measured 
attenuation  is  feasible.  We  would  not  anticipate  quantitative  agree¬ 
ment  because  of  the  air-earth  interface.  9 


*The  dielectric  constant  and  conductivity  of  the  medium  are  assumed  to 
be  constant  with  frequency  and  the  shape  of  the  curves  vs  frequency 
will  change  with  changes  in  the  assumed  load  impedance  of  the  receiving 
antenna.  The  results  in  Fig.  6  assume  a  constant  load  of  50  ohms.  For 
a  change  of  path  length  ah,  the  total  attenuation  for  ah  feet  is 
V, 


20  logio(y^)  where  V*2  and  V£l  are  the  voltages  for  paths 
Hi 

respectively  and  ah  =  h-j -H2- 
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Fig.  4.  Received  voltage  for  transmission  between  two  wire  antennas 
in  a  homogeneous  conducting  medium,  1  volt  input. 


CO 

ru 


(U 

4-> 


ro 


<U 

1.  • 

S  3 

CL 

o  c 

2  •<- 

+J 


<u  o 

<u  > 

+J  I— 
OJ 

_Q  * 

E 


O  *r- 

•I-  x> 

co  a; 

CO  E 


E  CT) 

co  c 

C  "i- 
ro  4-> 

l-  u 


J-  c 

°  8 


<U  co 
CTi  3 

to  O 

+o  <U 


O  <V 
>  CT> 
O 


■§ 


>  4= 


<U  10 
O 

(U  c 
cc  •(- 


LD 


CO 


! 

i 


i 

l 

17 

I 

I 


/ 


It  should  also  be  noted  that  the  wire  antenna  program  was  used 
to  predict  the  input  impedance  of  the  buried  antenna.  Direct  re¬ 
flection  measurements  on  the  buried  antenna  demonstrated  that  the 
predicted  value  was  within  ±10%  of  the  measured  value. 

Clearly,  the  comprehensive  analysis  and  compute  program  can 
be  used  to  generate  a  variety  of  probe  design  curve:  for  given 
target  ranges,  media  and  frequency  ranges.  As  noted  earlier,  we 
have  deferred  such  an  extensive  computational  program  in  favor  of  an 
experimental  program. 

Certain  studies  initiated  during  the  first  half  of  the  contract 
and  reported  in  the  Semiannual  Technical  Report  (see  Appendix)  were 
not  continued  during  the  second  phase.  The  grounded  monopole  antennu 
used  as  a  transmitter  to  make  certain  propagation  tests  was  discarded 
for  tunnel  hazard  detection  applications  for  two  basic  reasons.  First, 
this  unbalanced  antenna  as  a  receiver  has  very  severe  noise  problems. 
Second,  it  was  found  that  the  fields  produced  by  the  grounded  mono¬ 
pole  were  confined  to  very  narrow  depths  (essentially  the  depth  of 
the  mcnopole)  and  have  a  null  along  the  axis  of  the  monopole.  While 
the  grounded  mor  *,ole  is  not  appropriate  for  tunnel  problems,  the 
propagation  modt  uHng  the  monopole  as  a  transmitter  and  a  balanced 
dipole  as  a  rece  ?t.-r  would  appear  to  have  possible  application  for 
constitutive  parameter  measurements  at  shallow  depths.  An  analysis 
of  a  line  source  above  a  half-space  was  also  not  continued.  At  the 
suggestion  of  the  Technical  Project  Officer,  our  efforts  were  focussed 
on  the  experimental  program.  The  limited  budget  made  it  impossible  to 
complete  both  theoretical  and  experimental  programs.  Ultimately,  it 
will  be  extremely  important  to  evaluate  the  fields  of  various  sources 
in  functional  form  in  order  that  the  physical  mechanises  can  be  under¬ 
stood.  Numerical  analyses  such  as  that  used  for  the  wire  antenna  do 
not  yield  thit  information.  The  results  reported  in  the  Semiannual 
Report  (see  Appendix)  were  for  the  transient  fields  produced  by  a 
line  source  above  a  conducting  half-space  at  locations  either  directly 
below  the  source  or  on  the  surface  of  the  half-space.  The  dielectric 
constant  and  conductivity  of  the  conducting  medium  are  assumed  to  be 
constant  with  frequency  and  both  short  and  long  time  results  are  given. 
That  is,  in  the  limit  of  high  and  low  frequencies  the  wave  number  of 

the  conducting  medium  (4 (a+joje ) )  goes  to  jfuAje  and  fjwuo  respectively. 
Analytical  expressions  for  the  transient  fields  for  these  cases  are 
given.  It  may  be  possible  to  extend  these  results  to  an  arbitrary 
observation  point  within  the  medium  using  numerical  integrations.  It 
does  not  appear  feasible  to  consider  a  frequency-dependent  dielectric 
constant  or  conductivity  for  the  medium  even  if  a  mathematical  model 
for  the  parameters  were  available  -  the  analytical  forms  of  the  inverse 
transformations  would  become  hopeless  and  render  this  approach  invalid. 
This  does  not  negate  the  result  but  does  dictate  caution  in  their 
interpretation. 
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V.  ELECTROMAGNETIC  PULSE  SOUNDING  SYSTEM. 


A  block  diagram  of  the  electromagnetic  pulse  sounding  system 
with  a  direct  interfacing  to  the  computer  is  shown  in  Fig.  8  together 
with  one  version  of  the  electromagnetic  probe.  The  required  instru¬ 
mentation  modifications  for  remote  site  measurements,  i.e.,  divorcing 
the  system  from  a  direct  interfacing  with  the  computer,  will  be  dis¬ 
cussed  later.  The  probe  shown  in  Fig.  8  consists  of  two  grounded, 
orthogonally  oriented  dipoles,  where  each  dipole  arm  is  a  U  arrange¬ 
ment  of  two  wires.  The  probe  dimensions  shown  are  appropriate  for 
interrogation  of  targets  to  depths  of  30  to  40  feet  in  the  overburden 
using  the  near-zone  fields  of  the  probe.  For  interrogation  to  the 
sane  depths  in  a  different  medium,  the  overall  dimensions  of  the  probe 
will  remain  constant  but  the  shape  of  the  dipole  arms  would  be 
changed.  If  the  medium  is  such  that  the  use  of  ground  rods  is  not 
practical,  then  an  alternative  arrangement  using  metal  plates  to 
ohtain  a  capacitive  coupling  could  be  used.  The  purpose  of  the 
ground  rods  or  plates  is  tu  match  the  low  frequency  energy  in  the 
interrogating  pulse  into  the  ground.  In  the  next  section  the  tuning 
achieved  by  both  changes  in  the  shape  of  the  dipole  arm  and  ad¬ 
justments  of  the  ground  rod  depths  is  demonstrated. 

The  pulse  sounding  system  can  be  operated  in  two  modes;  a 
direct  mode  where  one  dipole  is  used  for  transmission  and  reception 
and  an  orthogonal  mode  where  transmission  is  on  one  dipole  anl  reception 
on  the  orthogonal  dipole.  In  the  orthogonal  mode,  the  probe  is  "blind" 
to  homogeneous  targets  including  the  air-medium  interface  and  "sees" 
only  inhomogeneous  targets.  The  direct  mode  "sees"  both  homogeneous 
and  inhomogeneous  targets.  It  should  not  be  inferred  that  the 
orthogonal  mode  will  not  "see"  planar-type  contrasts  of  infinite 
extent.  If  the  contrast  plane  has  some  tilt  with  respect  to  the  plane 
of  the  probe,  the  orthogonal  mode  will  indeed  respond.  Thus  the  only 
planar  discontinuity  for  which  the  probe  is  "blind"  is  one  parallel 
to  the  probe  plane  -  a  seemingly  remote  possibility  in  tunneling.  In 
operation  as  a  hazard  detection  device  in  tunneling  operations  it  is 
intended  that  for  each  location  of  the  probe  on  the  tunnel  face,  three 
measurements  will  be  made  -  two  direct  modes  using  the  two  orthogonal 
dipoles  individually  and  one  orthogonal  mode.  This  should  provide 
excellent  coverage  without  the  necessity  of  moving  the  probe.  During 
this  contract  period,  effort  has  been  primarily  concentrated  on 
orthogonal  mode  operation.  This  mode  affords  excellent  transmitter- 
receiver  isolation  and  is  ideally  suited  for  the  type  of  targets  con¬ 
sidered.  Also,  true  direct  mode  measurements  using  high  voltage  pulse 
generators  will  require  some  circuitry  modifications  in  order  to  protect 
the  sampling  head  of  the  scope  from  the  high  voltage  pulse.  The 
necessary  circuitry  is  felt  to  be  within  the  state-of-the-art.  The 
baluns  shown  with  the  probe  in  Fig.  8  are  special  broadband  baluns 
(Z-Line  50200E)  which  transform  an  unbalanced  coaxial  line  with  an 
impedance  of  50  ohms  to  a  balanced  transmission  line  with  an  impedance 
of  200  ohms.  A  variety  of  such  baluns  are  available  for  low  power 
applications  as  "off  fhe  shelf"  items.  However,  for  high  power  appli¬ 
cations  the  balun  design  to  achieve  broadbandedness  is  critical  and 
must  be  individually  tailored  for  specific  pulse  generators. 
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(b)  BLOCK  DIAGRAM  —  ORTHOGONAL  MODE 


Fig.  8.  (a)  Electromagnetic  probe. 

(b)  Block  diagram  of  pulse  sounding  system. 
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The  electromagnetic  pulse  sounding  system  currently  in  use 
employs  a  pulse  generator  which  produces  a  periodic  train  of  video¬ 
type  pulses.  The  pulses  are  coupled  into  the  medium  by  the  transmit 
arm  of  the  broadband  probe.  Periodically,  a  wave  is  launched  on  the 
transmit  dipole  at  the  feed  point  which  penetrates  into  the  medium. 

As  the  wave  propagates  along  the  dipole,  its  fringe  fields  penetrate 
deeper  into  the  medium.  The  dipole  is  not  currently  being  employed 
as  a  radiating  structure  in  the  normal  sense  since  it  is  the  near 
fields  which  are  being  used.  For  this  reason  the  dipole  is  referred 
to  as  a  probe  rather  than  as  an  antenna  structure.  Reflected  signals 
from  the  medium  are  observed  on  a  sampling  scope  acting  as  a  receiver. 
The  measurements  are  automated,  i.e.,  the  computer  controls  the  scope* 
sweep  and  records  a  response  waveform  which  is  an  average  of  some  200 
reflected  waveforms.  The  purpose  of  averaging  is  to  minimize  noise 
problems.  Once  the  averaged  waveform  is  recorded  it  is  available  for 
further  processing  either  in  real  time  or  at  some  later  time.  When 
the  optimum  processing  schemes  are  fixed,  the  computer  functions  can 
be  replaced  by  simple  circuitry,  but  at  this  research  stage  the  com¬ 
puter  offers  an  irreplacable  versitility. 


The  system  modifications  necessary  for  remote  measurements  are 
shown  in  Fig.  9.  An  external  ramp  generator  controls  the  horizontal 
axis  of  the  oscilloscope  display  providing  a  sweep  rate  independent 
of  the  sampling  time  base.  The  ramp  generator  also  provides  a  gate 
signal  which  goes  positive  when  the  sweep  begins  and  negative  when 
the  sweep  is  finished.  The  gate  and  ramp  signals  are  recorded  in 
parallel  with  the  oscilloscope  waveform  on  three  channels  of  a 
multichannel  FM  tape  recorder.  In  playback,  the  three  channels  recorded 
are  connected  to  three  analog/digital  converter  inputs  for  conversion 
of  the  oscilloscope  waveforms  into  digital  information. 


The  conversion  technique  is  Illustrated  in  Fig.  10.  The  ramp 
voltage  at  any  point  is  proportional  to  the  deflection  across  the 
screen  of  the  oscilloscope  waveform  that  has  been  recorded.  Thus  n 
evenly  spaced  samples  may  be  taken  by  taking  a  sample  at  the  start 
of  the  sweep  and  then  monitoring  the  ramp  and  taking  another  sample 
each  time  the  ramp  voltage  increases  by  aV  =  (Vo-Vi )/(n-l )  until  n 
samples  have  been  taken.  The  gate  voltage  is  used  to  provide  a 
precise  time  mark  for  the  first  sample.  From  then  on  the  sampling 

foHnfn*r^<  e5?n?in?d  by  the  ramp  volta9e-  This  ^formation  is  then 
ted  into  the  digital  computer  and  processed  in  the  same  manner  as 
data  from  the  local  measurements. 


A.  Probe  Configurations  and  Tests 

A  sketch  of  one  version  of  the  probe  structure  suitable  for 
interrogatKin  of  targets  at  depths  of  30  to  40  feet  was  shown  in 
Mg.  8.  In  Figs.  11  and  12  the  tuning  mechanisms  afforded  by, 

rnHceat1V^iii*  I he+ S5ape-r?f  thp  Jipole  arms  a"d  depth  of  the  ground 
rods  are  illustrated.  The  results  shown  in  Figs.  11  and  12  are  direct 

reflectometry  measurements  on  one  dipole  over  the  local  overburden 
h!  1J?ldent  Puls®  is  a  step  function  with  a  2  volt  peak  and  a 
.  "c  The  advanta9es  of  a  d>pole  arm  composed  of  two 

wires  in  a  U  or  V  arrangement  as  compared  to  an  arm  made  up  of 

*  wire  are  clearly  evident  in  Fig.  11  where  reflection  wave¬ 
forms  for  these  dipoles  are  compared  with  those  with  an  open  circuit, 
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REMOTE  MEASUREMENT  SYSTEM  -  RECORDING 


Fig.  9.  Remote  measurement  system. 

(a)  recording 

(b)  playback. 
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short  circuit  and  matched  load  at  the  probe  feed  point.  In  Fig.  12, 
the  tine  scale  is  compressed  in  order  to  illustrate  the  tuning 
capabilities  of  the  ground  rods.  By  adjusting  the  ground  rod  depth, 
the  later  portion  of  the  waveform  can  be  matched  to  the  waveform 
obtained  when  a  matched  load  is  placed  at  the  end  of  the  feed  cable. 

Note  that  the  waveforms  shown  in  Figs.  11  and  12  are  ample  evidence 
that  low  frequency  energy  can  be  effectively  coupled  into  the  medium. 
When  the  probe  configuration  of  Fig.  8  was  tested  in  the  orthogonal 
mode  it  was  found  that  the  system  was  extremely  insensitive  to  surface 
obstructions.  A  station  wagon  driven  between  the  probe  arms  had  no 
effect  on  the  response  waveform. 

The  probe  configuration  in  Fig.  8  was  designed  for  use  with  an 
HP  214-A  pulse  generator  purchased  for  use  on  this  contract.  The 
characteristics  of  the  HP  214-A  pulser  (50  volt  peak,  45  ns  base,  10 
KHz  repetition  rate)  are  appropriate  for  interrogation  of  subsurface 
targets  with  a  maximum  linear  dimension  of  less  than  10  times  the 
repetition  wavelength  in  the  medium.  Two  other  pulse  generators;  a 
laboratory  built  unit  (5  volt  peak,  3  ns  base,  1-2  MHz  repetion  rate) 
and  an  Ikor  generator  (1000  volt  peak,  250  ps  base,  250  Hz  repetition 
rate)  are  presently  available  at  this  laboratory.  The  Ikor  unit  be¬ 
came  available  less  than  one  month  ago. 

After  completion  of  a  series  of  design  measurements  similar  to 
those  shown  in  Figs.  11  and  12,  a  series  of  direct  and  orchoqonal 
mode  measurements  over  the  local  overburden  were  initiated.  It  was 
found  that  when  the  Hewlett  Packard  pulser  (HP214-A)  was  operated  at 
peak  power  (50  volt  peak  pulse),  the  balun  (50200E)  on  the  transmit 
probe  saturated.  Saturation  occurred,  in  fact,  for  pulse  peaks  less 
than  10  volts.  Saturation  of  the  transmit  balun  results  in  the  ap¬ 
pearance  of  false  targets  in  the  response  waveform  whose  position 
is  dependent  on  the  incident  pulse  magnitude  and  duration.  The  balun 
problem  does  not  in  itself  represent  a  basic  limitation  since  baluns 
with  a  higher  power  handling  capability  can  be  built.  At  the  same 
time,  one  should  not  underestimate  the  design  difficulties  inherent 
in  building  a  balun  to  handle  the  HP  214-A  pulse  generator  operating 
at  peak  power  and  to  achieve  a  bandwidth  from  10  KHz  to  60  MHz.  We  seek 
to  achieve  reasonable  pulse  fidelity  so  that  simply  being  able  to 
force  some  power  through  the  balun  over  this  frequency  range  is  not 
sufficient.  It  was  decided,  primarily  because  balun  winding  remains 
essentially  an  art  rather  than  a  science,  to  order  appropriate  baluns 
from  a  commercial  specialist  (Z-line  Corporation)  and  circumvent  the 
immediate  difficulty  either  by  eliminating  the  need  for  the  balun  or 
by  using  a  lower  power  pulse  generator.  The  decision  not  to  attempt 
to  wind  our  own  baluns  appeared  to  be  a  wise  one  since  even  for  those 
with  expertise  in  this  area,  a  promised  delivery  date  of  2  months 
extended  to  almost  5  months  indicating  some  design  difficulties.  We 
now  have  baluns  which  meet  the  specifications  of  the  HP214-A  pulse 
generator.  Unfortunately,  the  baluns  were  delivered  too  late  to 
attempt  a  measurement  sequence  on  subsurface  targets  during  this 
contract  period. 

Attempts  to  eliminate  the  need  for  a  balun  were  unsuccessful. 

In  each  scheme  tested,  it  was  found  that  when  the  baluns  were 
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removed  currents  on  the  feed  cables  produced  a  sensitivity  to  feed 
cable  position  which  was  intolerable.  Consequently,  a  series  of 
measurements  using  a  smaller  version  of  the  probe  and  the  laboratory- 
built  pulse  generator  were  initiated.  With  a  1  to  2  MHz  repetition 
rate  and  a  5  volt  peak,  this  pulse  generator  is  not  appropriate 
for  deep  penetration  of  the  medium.  However,  available  50200E 
baluns  pass  this  pulse  with  good  fidelity  and  thus  the  feasibility 
of  the  pulse  sounding  scheme  could  be  tested  on  shallow  targets, 
depths  to  5  feet  or  so  In  the  overburden,  without  encountering  balun 
saturation  problems.  Accordingly,  a  series  of  measurements  on  con¬ 
ducting  and  dielectric  cylindrical  targets  burled  to  depths  of  5 
feet  In  the  local  overburden  were  made.  Results  of  these  measure¬ 
ments  are  given  in  Section  VI.  The  probe  structure  used  for  these 
measurements  is  shown  in  Fig.  13.  Note  that  a  section  of  300  ohm 
twin  lead  has  been  Inserted  between  the  balun  output  and  the  feed 
point  of  the  dipoles.  The  purpose  of  the  twin  lead  is  to  clean-up, 
as  far  as  possible,  the  range  window  for  targets  at  shallow  depths 
In  the  overburden.  A  direct  or  reflection  mode  measurement  on  one 
dipole  of  the  probe  in  Fig.  13  Is  shown  In  Fig.  14.  Note  that  the 
balun  reflection  is  primarily  high  frequency  and  could  be  anticipated 
since  the  high  frequency  end  of  the  50200E  balun  is  around  60  MHz. 

A  second  effect  was  noted  from  measurements  similar  to  the  one  In 
Fig.  14.  It  was  found  that  changing  the  ground  rod  depth  had  little 
effect  on  the  response  waveform.  Comparing  this  result  with  those 
obtained  for  the  larger  prob  Indicates  that  tuning  via  ground  rod 
depth  adjustment  Is  primarily  effective  for  frequencies  below  1  MHz, 
at  least  when  the  probe  is  located  on  the  overburden.  On  a  less 
lossy  medium,  sensitivity  to  the  ground  rod  depth  should  extend 
to  higher  frequencies. 

Early  In  the  contract  period  it  was  decided  that  some  remote 
site  measurements  should  be  made  at  a  nearby  limestone  quarry. 

The  limestone  medium  (cr  %  8.0,  o  <  10-5  nhos/meter)  was  the  most 
practical  approximation  to  a  hard  rock  medium  (e.g.,  granite, 
cr  1.  7.0,  o  <  lO-o)  within  a  reasonable  distance  from  the  laboratory. 
Accordingly  a  mobile  rig  was  outfitted  (power  to  run  the  pulse  gen¬ 
erators  and  scopes  was  not  available  at  the  site)  and  Instrumentation 
to  divorce  the  system  from  the  computer  devised.  Certain  man-made 
and  geological  anomalies  at  the  quarry  appropriate  for  Interrogation 
with  the  large  probe  and  HP214-A  pulse  generator  had  been  selected. 
When  It  became  evident  that  the  high  power  baluns  would  not  be 
available  during  this  contract  period,  It  was  decided  to  take  a 
series  of  measurements  at  the  quarry  using  the  smaller  probe  and 
pulse  generator.  It  was  known  In  advance  that  certain  of  the  targets 
selected  were  not  appropriate  using  this  pulse  sounding  system. 

The  results  of  these  measurements  are  given  In  Section  VII. 

Using  flexible  wires  to  construct  the  arms  of  the  probe  has 
the  advantage  that  the  wires  can  be  made  to  conform  to  the  Individual 
undulations  of  the  medium.  A  major  disadvantage  is  that  aligning 
the  Individual  dipoles  to  achieve  precise  orthogonality  Is  often  a 
tedious  and  time-consuming  effort.  Alternatively,  constructing  the 
probe  arms  of  rigid  members  automatically  yields  orthogonality,  but 
now  the  height  of  the  feed  point  can  vary  and  the  probe  may  well  have 
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Fig.  13.  Small  electromagnetic  pulse  sounding  probe. 


Fig.  14.  Direction  reflection  mode  of  one  probe. 
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some  slope  with  respect  to  the  medium's  mean  surface.  Certain  tests 
have  been  made  of  small  probe  configurations  with  respect  to  optimum 
arm  configurations,  feed  point  heights  and  the  use  of  an  absorbing 
material  between  the  probe  and  the  overburden  (the  vegetation  above 
the  overburden  highly  attenuates  the  high  frequency  energy  in  the 
pulse).  The  results  of  these  tests  are  summarized  below. 

A  preliminary  study  of  probe  parameters,  including  effects  of 
varying  conductor  shape,  height  above  the  ground  surface,  and  placing 
absorptive  material  between  the  probe  and  ground  surface,  has  been 
made.  In  all  cases  an  orthogonal  transmit-receive  probe  pair  was 
used.  This  approach  seems  to  be  the  best  way  of  isolating  two  probes 
which  are  very  close  to  one  another. 

When  a  probe  pair  lies  on  a  rough  surface,  it  is  known  that 
isolation  is  degraded  due  to  the  asymmetry  caused  by  the  uneven 
ground.  The  frequency  content  of  this  coupling  is  inversely  pro- 
po 'tional  to  the  average  sizes  of  the  ground  irregularities,  but  in 
general  lies  above  40  MHz.  For  an  application  concerned  with  distant 
targets,  where  return  with  frequency  content  below  about  40  MHz  is 
expected,  the  ground  irregularity  problem  is  of  little  concern,  since 
it  can  be  virtually  eliminated  by  filtering.  It  was  found  that 
raising  the  probes  2  to  3  inches  above  the  ground  and  using  ab¬ 
sorptive  material  between  the  probe  and  ground,  both  improved  the 
target  return  to  clutter  ratio  for  shallower  targets  where  hiqher  frpouencv 
returns  were  expected.  J 

The  effects  of  antenna  geometry  and  conductor  shape  are  seen  at 
al  frequencies.  For  the  "bow-tie"  (Fig.  15)  configuration  used  in 
this  study,  the  ends  of  the  structure  have  been  identified  as  the 
major  source  of  coupling  problems.  Asymmetries  in  layout  or  ground 
parameters  at  the  ends  result  in  a  degraded  target  return  to  clutter 
ratio.  Initial  tests  have  shown  that  the  elimination  of  multiple  re¬ 
flections  on  the  probe  structure,  and  a  "smoother"  shape  in  general  result 
in  an  improved  performance.  Further  tests  using  models  appropriate  for 
distant  targets  are  planned. 

B.  Processing  Pulse  Sounding  Measurements 

Simplified  block  diagrams  for  the  electromagnetic  pulse  sounding 
probe  in  the  orthogonal  mode  are  shown  in  Fig.  16  for  no-target  (Fig. 

16a)  and  target  (Fig.  16b)  configurations,  where  the  block  T  indicates 
the  presence  of  the  target.  Capitalized  letters  correspond  to  fre¬ 
quency-dependent  quantities  and  the  lower  case  counterpart  to  a  time- 
dependent  waveform.  Each  capital  and  lower  case  symbol  are  a  transform 
pair.  In  terms  of  the  complex  variable  S  =  a  +  ju 

ENT(S)  *  Ct(S)  Pt(S)  G(S)  Pr(S)  Cr(S)  e1(S> 

(1) 

E?(S)  =  Et(S)  Pt(S)  6(S)  T(S)  Pr(S)  Cr(S)  Ei(S) , 
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Fig.  15.  Sketch  of  experiment  probe  geometry. 


(a) 


(b) 


Fig.  16.  Block  diagram  or  pulse  sounding  probe 
in  orthogonal  mode. 

(a)  no-target 

(b)  target. 
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where  C,  P,  G  and  T  are  the  transfer  functions  of  the  cable,  probe 
ground  and  target  respectively  and  the  subscripts  t  and  r  denote 
transmit  and  receive.  The  cable  transfer  functions  include  the 
baluns  and  short  lengths  of  twin-lead  line.  Making  the  reasonable 
assumptions  that  the  transmit  and  receive  properties  are  identical 
Eq.  (1)  simplifies  to  ’ 

=  C(S)  P(S)2  G(S)  E1  (S) 

(2) 

E?(S)  =  C(S)2  P(S)2  G(S)  T(S)  Ei(S)  . 

By  shorting  the  cables  together  including  the  baluns  and  twin-lead 
line,  the  normalizing  response 

(3)  E{j(S)  =  C(S)2  E1' (S) 
can  be  recorded.  Ideally  then 

E°(S)  -  EfjT(S) 

(4)  — 7o/0 - =  P(S)2  G(S)  [T(S)  -  1]  . 

The  normalizing  response  e2(S)  is  stable  and  repeatable.  However, 
the  no-target  measurement  is  not  precisely  stable  if  the  probe  is 
moved  to  various  ground  locations.  This  is  due  to  localized 
properties  of  the  ground  such  as  roughness,  inhomogeneities, 
moisture  content,  vegetative  cover,  etc.  which  changes  G(S).  Some 
is  al so  due  to  changes  in  P(S)  which  is  somewhat  sensitive 
to  the  height  of  the  probe  feed  points  above  the  ground.  We  choose 
to  ignore  the  probe  effect  since  most  of  this  occurs  in  the  direct 
coupling  between  the  transmit  end  recei vie  probes  and,  as  will  be 
seen  later,  can  be  gated  out.  At  the  jth  location  then  the  no¬ 
target  response  is  more  properly 


(5)  Ej}J(S)  =  C(S)2  P( S ) 2  E1  (S)  CGq(S)  +  ^'(S)]  . 


In  practice,  some  M  no-target  measurements  are  made  and  an 
average 

_ 

(6)  ent(s)  *  M  .1  eSt(s)  -  CIS)2  P(S)2  e'(S)  CG0(S)  +  i  [  NJ(S)] 

J-l  n  j=l 

is  obtained.  In  a  similar  fashion  an  average  target  waveform 
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(7) 


same  Pr'>be-targetrgec»netrynorabynsimSivXDutt^d  keep1n9  the 

times  in  the  same  location  Nn+a  JhPly«.[!utt1n9  the  Probe  down  M 
has  been  writtln  as  1,016  that  the  9round  transfer  function 


(8)  G(S)  =  G0(S)  +  nj(s), 

for  the  no-target  case  and  as 


(9)  G(S)  =  Gq(S)  +  NJ(S) 


for  the  target  case.  The  reason 
different  in  that  for  the  target 
The  normalized  difference  of  the 


being  that  the  sampling  is  actually 
case  the  variations  are  restricted 
averages  is 


(10) 


E?(S)  - 


eJt(S) 


e£(S) 


•  p(sy 


G0(S) 


CT(S)  -  1> 


If  we  make  the  assumption  that  N  =  N  then 

01) 

E?(s)  -  eJt(s) 

F0,F.  ■  (T(S)-DP(S)2 

Let 

1  M  1 

M  X  N^S)  -  N(S)  , 

J  ™  1 

(12) 

then 

(13) 

E?(S)  -  E°t(S) 

FoK ,  -  (t(s)-dp(s)2 

or5' 

re0{s) 


1  M  • 

+  »  X  ^(sn 
* 


CGq(S)  +  N(S)]  . 
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Interesting  comparisons  to  Eg.  (13)  are 


(14) 

and 

(15) 


E$(S)  -  E°t(S) 


NT 


E„(S) 


E°J(S)  -  E“(S) 


■  P(S)2  Itt^S)  -  N(S)-| 


E“(S) 


P(S)2  T(S)  [N%)  -  N(S)3  . 


which  essentially  show  the  variance  of  the  no-target  and  tarqet 
measurements  respectively. 

Actually,  the  differencing  and  averaging  are  done  in  the  time 
domain 


The  fast  Fourier  transform  of  the  difference  waveform  is  then  normalized 
by  the  pulse  and  cable  spectrum.  It  has  been  demonstrated  that  the 
transform  of  the  difference  waveform  is  the  difference  of  the  trans¬ 
formed  waveform,  thus  no  aliasing  errors  are  present.  Usually,  the 
difference  waveform,  <5(t),  is  gated  before  being  transformed.  That  is, 
the  difference  waveform,  6(t),  is  multiplied  by  a  gate  waveform 


(17)  g(t)  *  a(t-T] )  Qi(t-T.| )  -  u(t-T2)3  + 

1  Cu(t-T2)  -  u ( t-T 3 ) D  + 

□  -  a(t-T3)]  Cu(t-T3)  -  u(t-T4)3  , 

where  a,  Tlt  T2,  etc.  are  real  constants.  The  time  durations  T?-Ti 
and  T4-T3  are  equal,  and  can  be  zero.  The  purpose  of  the  gate  is  to 
remove  the  leading  portion  of  the  difference  waveform  which  has  rather 
high  frequency  content  and  can  be  attributed  to  direct  coupling  be¬ 
tween  the  transmit  and  receive  probes.  It  can  also  be  used  to  isolate 
the  scattered  fields  from  different  reflection  mechanisms.  Examples 
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of  the  ungated  and  gated  difference  waveforms  will  be  shown  liter. 

A  trapazoidal  gate  is  used  in  order  to  avoid  placing  artifically 
generated  step  discontinuities  in  the  gated  waveform. 

Processing  the  data  then  involves  the  following  steps: 

1)  averaging  the  target  and  no-target  time  domain  waveforms. 

2)  Obtaining  the  time  domain  difference  waveform,  i.e.,  the 
target  average  minus  the  no- target  average. 

3)  Gating  the  difference  waveform  to  remove  the  direct  couplinq 
between  transmit  and  receive  probes. 

4)  Transformations  of  the  gated  difference  waveform  to  the 
frequency  domain,  i.e.,  a  fast  Fourier  transform. 

5)  Normalization  of  the  transformed  difference  by  the  pulse- 
cable  transfer  function. 

6)  Inverse  transformation  of  the  normalized  spectrum  in  5) 
back  to  the  time  domain  -  inverse  fast  Fourier  transform. 

Four  signatures  of  the  target  are  available  for  study;  the 
amplitude  and  phase  characteristics  in  the  frequency  domain  (5  above), 
the  time  domain  difference  waveform  before  normalization  (3  above) 
and  the  time  domain  difference  waveform  after  normalization  (6  above). 
From  a  comparison  of  these  signatures  for  various  targets,  one  can 
begin  to  consider  the  separation  of  targets  and  the  recognition  of 
a  given  target.  We  have  confined  our  results  for  the  present  to 
examinations  of  the  normalized  amplitude  spectrum  and  the  gated 
difference  waveform  before  normalization.  The  phase  spectrum  does 
have  information  content  but  it  is  generally  more  difficult  to  utilize. 
Also,  we  will  generally  need  more  power  and  lower  frequencies  in  the 
interrogating  pulse  spectrum  in  order  to  obtain  reasonable  phase 
spectrums.  With  more  advanced  processing  schemes  we  would  begin  to 
comtine  the  amplitude  and  phase  data  at  several  frequencies  to 
synthetically  produce  waveforms  uniquely  characteristic  of  the 
target  and  possibly  intimately  related  to  the  physical  properties  of 
the  target.  As  was  demonstrated  earlier,^  obtaining  waveforms 
predictably  dependent  upon  the  target's  physical  properties  may  not  be 
possible  in  a  highly  dispersive  medium  but  this  does  not  negate  ob¬ 
taining  a  characteristic  waveform. 

In  the  next  two  sections  of  this  report,  certain  of  the  wave¬ 
forms  and  spectra  d-scussed  above  for  cylindrical  targets  in  the 
overburden  and  for  targets  in  a  limestone  medium  are  presented. 
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C .  Manipulation  of  Time  Domain  Ref 1 ectometry  Data 

It  was  noted  earlier  (Section  IVB)  that  an  extremely  versatile 
program  for  the  automated  measurement  and  processing  of  pulse 
sounding  data  was  used  with  a  directly  interfaced  instrumentation 
computer.  This  software  has  been  developed  on  a  number  of  programs 
including  this  contract.  The  main  features  of  the  program  are 
described  below. 

Program  for  Manipulation  of  Time  Domain  Data 

This  program  provides  a  means  of  recording  time  domain  data 
from  a  sampling  oscilloscope,  drawing  it  on  the  plotter,  storing  it 
on  the  disc  and  performing  various  manipulations  including  Fourier 
transforms. 

Data  files  consist  of  a  single  disc  track  and  must  be  created 
with  the  DAT  command  prior  to  program  execution.  They  consist  of 
256  data  words  followed  by  10  file  description  words  and  finally  a 
code  word  which  is  equal  to  integer  1  for  a  frequency  domain  file 
and  zero  for  a  time  domain  file. 

The  program  actually  consists  of  three  complete  programs  which 
are  interchanged  by  the  core  swapping  subroutine  SWAP.  To  load  the 
program: 

First  load  TDRSS  in  the  regular  way  and  execute. 

Then  load  TDRS  in  the  regular  way  and  execute. 

Finally  load  TDRM  in  the  regular  way  and  execute. 

The  program  has  its  own  operating  system  and  when  TDRM  is 
executed  it  will  respond  with  a  #  sign.  It  will  then  accept  any  of 
the  following  commands  (only  the  first  three  characters  are  sig¬ 
nificant). 

CALibrate  -  The  sampling  'scope  is  set  to  the  uppe’-  scale  limit 
and  a  carriage  return  is  typed.  Repeat  for  lower 
scale  limit.  Computer  gives  instructions.  These 
points  will  be  ±4  inches  on  the  plotter. 

RECord  -  256  samples  are  recorded  from  the  sampling  'scope 
(average  of  200  sweeps)  and  converted  to  a  ±4  scale 
in  terms  of  the  calibration  points. 

WRIte  file-  Writes  the  256  data  points  currently  in  core  on  the 
disc.  Program  responds  with  NAME:  The  file  name  is 
then  typed  in.  Program  then  prints  DESCRIPTION.  Up 
to  30  characters  of  description  may  then  be  typed  in 
which  are  written  on  the  file.  Program  follows  with 
INSTRUMENT  SETTINGS.  A  further  30  characters  may  be 
used  to  indicate  these.  The  file  is  written  after 
the  final  carriage  return. 
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READ  file  -  NAME  is  requested  and  the  previously  written  file  read 
back  into  core. 


DEScription-Prints  out  the  30  character  description  and  the  30 
characters  indicating  instrument  settinqs  of  the 
most  recently  read  or  written  file. 


WHAT  is?  - 


PLOT 


XSHift 


YSHift 


GRId 

ADVance 

REVerse 

FFT 


IFFt 

Filter 


Program  requests  NAME.  Thr  30  character  description 
of  the  named  filejs  then  printed  out  without  affectina 
thema^core  buffer.  Used  to  check  wHFt'h'er  an  old — S 
tile  is  wanted  before  writing  over  it. 

■  nira^  the  rurrfnt  contents  of  the  core  buffer  on  the 
plotter.  In  the  time  domain  the  scale  is  simply  +4 

I"  the  !re?uency  domain  only  the  amplitude  of 
a  t  to -SoTb  s«le  f’rSt  1M  ha™°n1CS  15  P'0tted  on 

'  KL’S  r^uest!d-  The  data  PoiP‘5  are  then  shifted 
in  core  by  the  indicated  amount.  Points  are  shifted 

Points  nn^  12CaJi?U;  there  is  no  interpolation. 

'?"ts  ]°st  a*  one  end  of  the  trace  are  added  to  the 
other  end.  (Do  not  use  for  frequency  domain  data.) 

-  INCHES  in  DBS  are  requested  for  the  time  or  frequency 
domains  respectively.  The  data  in  core  is  then  shifted 
by  the  indicated  amount.  In  the  frequency  domain  only 

changed!  1s  shifted*  The  angles  remain  un- 

-  Draws  a  grid  on  the  plotter.  The  computer  requests 
the  number  of  divisions  required  in  each  direction. 

-  Moves  the  plotter  forward  one  page. 

-  Moves  the  plotter  back  one  page. 

'  corf  hiff, JaSt™°l,!;1?r  transfom  of  the  contents  of  the 
core  buffer.  The  data  output  is  stored  so  that  the  odd 

number  points  are  amplitudes  on  a  dB  scale  starting  with 

location  continuin9  to  the  127th  harmonic  in 

ocation  255.  The  even  nunber  points  are  phase  angles 

in  radians.  The  0  dB  level  corresponds  to  the  d.c  power 
doma?nent6d  ^  3  COnstant  4  inch  deflection  in  the’time 

Takes  the  inverse  Fourier  transform. 

nllt  J?lfI°CeST,a  tinie  d0main  wavefo™  through  a  band 
i4  *  -The  pr?gram  requests  the  hiqhest  and 
jowest  harmonics  to  be  saved  (0  to  127,  o' being  the 
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COMbine  -  Provides  for  adding  or  subtracting  files  in  any  manner. 

files  FILES  (NAME  AND  MULTIPLIER)  is  printed  followed  by  a 

col  an  (:).  The  file  name  is  typed  in  followed  by  a 
carriage  return.  An  X  is  then  printed.  Any  real 
multiplier  is  then  typed  in  followed  by  a  carriaqe 
return.  This  sequence  is  repeated  as  many  times  as 
desired.  A  carriage  return  in  response  to:  terminates 
the  sequence  leaving  the  weighted  sum  in  the  core  buffer 
(e.g.,  to  obtain  a  difference  curve  a  multiplier  of  1. 
could  be  used  for  the  first  file  and  -1.  for  the  second. 
Multipliers  of  5.  and  -5.  would  produce  a  magnified 
difference.  A  multiplier  of  0.25  for  each  of  four 
files  would  produce  an  average,  etc.)  Taking  a  difference 
of  two  frequency  domain  files  produces  a  ratio  since 
amplitudes  will  be  subtracted  in  dBs  and  phases  in 
radians.  (A  multiplier  used  here  would  amount  to 
raising  to  a  power).  This  permits  impulse  responses 
to  be  found  by  dividing  a  response  waveform  by  that 
of  an  incident  pulse;  also  reflection  coefficients  by 
dividing  a  response  by  that  of  a  short  circuit. 

PRInt  -  Prints  all  or  part  of  the  current  core  buffer  on  the 

teletype.  The  starting  element  number  will  be  requested. 
Printing  may  be  interupted  with  the  BREAK  button.  A 
control  D  will  then  return  to  the  #  sign.  The  scale  of 
the  print  out  is  in  inches  of  plotter  movement  in  the 
time  domain  and  dB  and  radians  in  the  frequency  domain 
(see  FFT  for  data  arrangement). 

LABel  -  Draws  the  current  DESCRIPTION,  INSTRUMENT  SETTINGS  and 
FILE  NAME  on  the  plotter.  LABEL  NUM:  will  be  requested. 
The  first  label  will  be  in  the  top  left  hand  corner  of 
the  page.  If  a  second,  third,  etc.  label  is  required  it 
will  be  moved  down  1"  each  time.  A  carriage  return  has 
the  same  effect  as  typing  number  1  for  the  label  number. 

GATE  -  Gates  out  a  portion  of  the  time  domain  waveform.  The 
GATE  OPEN  and  GATE  CLOSE  positions  In  Inches  are 
requested.  All  points  before  OPEN  and  after  CLOSE  are 
set  to  zero. 


VI.  PULSE  SOUNDING  DATA  FOR  CYLINDRICAL  TARGETS 
BURIED  IN  THE  OVERBURDEN 

Pulse  sounding  measurements  of  various  cylindrical  targets  buried 
in  the  overburden  close  to  the  ElectroScience  Laboratory  have  been 
made.  The  measurements  were  made  using  the  small  probe  shown  in  Fig. 

13  and  the  laboratory-built  (5  volt  peak,  3  ns  base,  1.5  MHz  repetition 
rate)  pulse  generator.  The  targets  consisted  of  2,  3  and  4  inch  diameter 
hollow  dielectric  (e  =  2.56)  pipes  at  a  depth  of  3  feet  in  the  over- 
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burden  and  4  inch  diameter  conducting  hollow  pipes,  1,  3,  5  and  10 
feet  long  at  a  depth  of  5  feet  in  the  overburden.  The  dielectric 
tarqets  were  some  20  feet  long  or  qreater.  Other  conducting  pipes 
in  the  same  sequence,  i.e.,  lengths  of  1 ,  3,  5  and  10  feet,  are  also 
buried  at  depths  of  1  foot  and  3  feet  and  have  been  measured,  but  the 
data  presented  here  will  be  confined  primarily  to  the  5  foot  depth. 

In  Figs.  17  and  18,  respectively,  the  interrogating  pulse  and 
its  spectrum  are  shown.  The  pulse  and  spectrum  shown  corresponds  to 
propagation  through  both  feed  cables  (roughly  150  feet  each),  both 
baluns  (50200E)  and  two  sections  of  300  ohm  twin  lead  shorted  together. 
The  balun  characteristics  shown  in  Figs.  17  and  18  will  not  change 
when  the  twin  lead  is  attached  to  the  probe.  The  balun  operation  is 
separated  in  time  from  the  probe  and  it  is  the  characteristic  im 
pedance  of  the  line,  not  the  total  impedance  seen  at  the  output  port 
of  the  balun,  which  is  important  for  the  short  time  response  of 
interest.  Clearly  the  pulse  shown  in  Fig.  17  is  not  the  waveform 
being  transmitted  into  the  ground  -  it  is  however  the  waveform  at  the 
feed  point  of  the  probe.  Ideally,  the  normalization  of  the  response 
waveforms  should  be  by  the  actual  transmitted  pulse,  but  this  waveform 
cannot  be  measured.  To  obtain  this  waveform  requires  an  absolute 
calibration  of  the  system  and  this  in  turn  requires  both  an  experimental 
measurement  of  a  known  target  geometry  and  an  analytical  solution  must 
be  obtained  using  the  near-zone  fields  of  the  probe.  Thus,  an  extention 
of  the  wire  antenna  analysis  is  needed  to  provide  a  complete  calibration 
of  the  system.  Note  that  the  spectrum  shown  in  Fig.  18  has  been 
normalized  by  setting  the  magnitude  of  the  harmonic  with  the  peak  arroli- 
tude  to  unity.  One  must  carefully  distinguish  between  this  type  of 
normalization  and  that  suggested  to  remove  the  Incident  pulse,  cables 
and  baluns  from  the  spectral  and  temporal  responses.  The  spectrum  in 
Fig.  18  corresponds  to  the  quantity  Ej^(S)  In  Eq.  (3),  where  it  is 
assumed  that  all  other  spectra  will  be  normalized  to  unity  first  in 
the  same  way. 


NANOSECONDS 


Fig.  17.  Interrogating  pulse  waveform. 
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^nra-cTf  in  Section  VB’  that  U  was  important  to  secure 

av  rages  of  both  the  no-target  and  target  response  waveforms  because 

of  possible  individual  variations.  In  Figs.  19a  and  19b  ™pectWelv 
dif^rreUnl^S«tionsthan9dn?1  TOde)  °f  l  no‘tar9et  confi guratioS  at  4  ' 

inch  a°;ash3o^ot 

orthogonal  data  reported  here  and  In  Section  VII  were  made  usi™  tuil 
con,  .gyration.  Shown  in  Fig.  19c  and  Fig.  19d  re“pect™elv  a^e  the 
average  waveforms  for  the  no-target  and  target  configurations  De- 

19ht|i?!d.fbKSei^  °"  c?"par,sons  of  Selected  pair  from  Figs  'l9a  and 

lutri::1:'  bui a  «"parison  »f 

e  ves  little  or  no  doubt  as  to  the  presense  of  a  target. 

SDectra  fir  Ihi’Jnn9  ?agef’  measured  and  Processed  wavefonns  and 
respectively^  "9  tar9ets  a^e  given  in  Figs.  20  to  28, 


1. 

2. 


8. 


10. 


inch  diameter  dielectric  pipe, 
inch  diameter  dielectric  pipe, 
inch  diameter  dielectric  pipe, 
inch  diameter  x  1  foot  length  conducting  pipe 

nch  d  ameter  x  3  foot  length  conducting  pipe 

nch  diameter  x  5  foot  length  conducting  pipe 

inch  diameter  x  10  foot  length  conducting  pipe 

dummy  trench  (5  foot  depth  x  10  foot  lenqth) 
i.e.,  a  trench  dug  and  then  refilled 
No-target  results. 


2 

3 

4 
4 
4 
4 
4 
a 


!anvffSnnsnan;espectt«Car“]hnownb!i”US  exceptions-  thp 


a. 

b. 

c. 

d. 

e. 


t*s  sverage  - 

the  unnormalized  amplitude  spectrum  of  the  gated 
difference  waveform. 

the  normalized  amplitude  specturm  of  the  gated 
difference  waveform. 

the  normalized  amplitude  spectrum  of  the  gated  difference 
an^baluns  removed aPeC‘rl,n  °f  th*  "“M“?  pP'”’ 

/  \  ^ev['ra^  general  points  should  be  noted.  The  difference  waveform 
?  Mh  9ate?  t0  bas1«"y  remove  the  direct  coupling  beSeen 

spectrin  ?d>  an  enhanc^  f" S  Mg^freg^cfS  r'e™  l*\?nT 
the  process  involves  division  by  small  ninbers  Pe.MsHrxiil  'x. 

300CMHa.  6  MCh  f,9ure  shou,d  be  ^regarded  above  roughi; 
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DIRECT  COUPLING  ACROSS  FEED  TERMINALS 


Fig.  19.  Orthogonal  mode  response  waveforms. 

(a)  no-target 

(b)  3  foot  cylinder  target 

(c)  no-target  average 

(d)  target  average. 


(b) 


Fig.  20.  Orthogonal  mode  response  2"  dielectric  pipe. 

(a)  difference  waveform 

(b)  gated  difference  waveform 

(c)  normalized  amplitude  spectrum 

(d)  normalized  amplitude  spectrum  with 
incident  pulse  removed. 
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Fig.  21.  Orthogonal  mode  response  3"  dielectric  pipe 

(a)  difference  waveform 

(b)  gated  difference  waveform 

(cj  unnormalized  amplitude  spectrum 

(d)  normalized  amplitude  spectrum 

(e)  normalized  amplitude  spectrum  with 
Incident  pulse  removed. 
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Fig.  23.  Orthogonal  mode  response  1'  cylinder. 

(a)  difference  waveform 

(b)  gated  difference  waveform 

(c)  normalized  amplitude  spectrum 

(d)  normalized  ampl  itude  spectrum  with 
incident  pulse  removed. 
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Fig.  24.  Orthogonal  mode  response  3'  cylinder. 

(a)  difference  waveform 

(b)  gated  difference  waveform 

(c)  unnormalized  amplitude  spectrum 

(d)  normalized  amplitude  spectrum 

(e)  normalized  amplitude  spectrum  with 
incident  pulse  removed. 
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Fig.  25.  Orthogonal  mode  response  5'  cylinder. 

(a)  difference  waveform 

(b)  gated  difference  waveform 

(c)  unnormalized  amplitude  spectrum. 

(d)  normalized  amplitude  spectrum 

(e)  normalized  amplitude  spectrum  with 
incident  pulse  removed. 
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Fig.  26.  Orthogonal  mode  response  10'  cylinder 
jaj  difference  waveform 
(b)  gated  difference  waveform 
(c|  unnormalized  amplitude  spectrum 
(aj  normalized  amplitude  spectrum 
(e;  normalized  amplitude  spectrum  with 
incident  pulse  removed. 


Fig.  27.  Orthogonal  mode  response  dummy  trench. 

(a)  difference  waveform 

(b)  gated  difference  waveform 

(c)  unnormalized  amplitude  spectrum 

(d)  normalized  amplitude  spectrum 

(e)  normalized  amplitude  spectrum  with 
incident  pulse  removed. 
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Fig.  28.  No-target  response,  orthogonal  mode. 

(a)  time  waveform  (average) 

(b)  gated  time  waveform 

(c)  amplitude  spectrum 

(d)  amplitude  spectrum  with  incident 
pulse  removed. 
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Consider  first  the  gated  difference  waveforms,  b,  for  each 
target.  It  is  clear  in  each  case  that  a  target  is  present.  This 
was  also  clear  when  the  average  target  response  was  compared  to 
the  average  no-target  response,  e.g.,  Figs.  19c  and  19d.  Note 
that  the  presence  of  the  dummy  trench,  Fig.  27,  is  also  evident. 

At  the  time  this  series  of  measurements  was  taken,  the  fciuried  targets 
had  been  in  place  for  approximately  3  months.  When  the  targets  were 
installed,  the  fill  was  compacted  in  each  case.  It  will  be  interesting 
to  observe,  on  a  continuing  basis,  how  long  the  dummy  trench  continues 
to  differ  significantly  from  the  no-target  response.  The  response 
of  the  dummy  trench  also  indicates  that  the  target  response  wave¬ 
forms  are  really  a  combination  of  the  target  and  trench  responses. 

Of  course,  no  such  dummy  trench  will  exist  in  advance  of  any  hard 
rock  tunneling  operation.  It  is  interesting  to  note  however  that 
perturbations  that  one  would  usually  assume  to  be  insignificant  are 
detectable. 

From  an  examination  of  the  gated  (b)  or  ungated  (a)  difference 
waveforms  the  following  conclusions  can  be  drawn; 

1)  the  class  of  dielectric  targets  can  be  separated  from  the 
class  of  conducting  targets,  and  both  can  be  separated 
from  the  dummy  trench  response. 

2)  From  the  difference  waveforms,  separation  of  the  conducting 
target  lengths  or  the  dielectric  target  diameters  does  not 
appear  feasible. 

3)  the  direct  coupling  between  transmit  and  receive  probes  (a) 
appears  significantly  different  for  the  dielectric  targets 
and  dummy  trench  than  it  does  for  the  conducting  targets. 

This  difference  is  the  result  of  the  way  the  individual 
measurements  to  secure  an  average  response  were  made.  Be¬ 
cause  the  dielectric  targets  were  extended,  the  individual 
measurements  were  made  by  moving  along  the  target  keeping 
the  same  probe-target  geometry.  For  the  conducting  targets, 
the  probe  was  simply  replaced  several  times  in  the  same 
location.  The  result  of  the  local  influence  of  the  grounds 
on  the  direct  coupling  is  particularly  noticeable  for  the 
perturbed  ground. 

4)  the  ungated  difference  waveforms  (a)  for  the  dielectric 
targets  endure  for  roughly  100  ns  and  the  waveforms  for 
the  conducting  targets  and  dummy  trench  for  approximately 
140  ns.  This  is  reasonable  considering  the  difference  in 
depth  of  the  targets.  It  also  indicates  that  in  the  case 
of  the  dummy  trench,  the  bottom  of  the  trench  is  seen; 
although  the  response  is  clearly  very  small. 
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A  casual  observation  of  the  normalized  spectra  of  the  gated  dif¬ 
ference  waveforms  with  the  incident  pulse  removed  would  indicate  that 
all  of  these  targets  appear  separable.  However,  as  will  be  discussed, 
the  influence  of  the  trench  is  the  major  source  of  the  differences, 
particularly  since  the  returns  from  the  trench  and  from  the  target 
are  not  completely  coherent. 

The  amplitude  spectra  for  the  gated  difference  waveforms,  with 
and  without  the  incident  pulse  removed,  show  much  the  same  type  of 
separations  and  distinctions  as  those  deduced  from  the  time  domain 
waveforms.  In  general,  however,  the  differences  for  target  classes 
are  much  less  distinct.  We  have,  however,  discarded  the  phase  spectrum. 
Clearly  both  trie  real,  time-dependent  response  waveform  and  its  complex, 
frequency-dependent  transform  contain  the  same  information.  Visually, 
in  these  cases,  the  time  waveforms  seem  more  simple  to  interpret. 

For  cylindrical  targets  of  length  L  we  would  anticipate  dipole- 
type  resonances  in  the  amplitude  spectra  at  the  following  frequencies 


f  =  hz  t 

"  2UFr 

which  yields  the  following  frequencies  in  MHz  for  the  conducting 
targets  assuming  c  =  10. 06 


n 

1  foot 

3  foot 

5  foot 

10  ft 

1 

158 

53 

32 

16 

2 

464 

159 

96 

48 

A  creeping  wave-type  resonance  {cylinder  radius  r) 


.  0.15 


for  r.  =10.0  yields  (for  a  4  inch  diameter) 
r 

f  =  286  MHz  and  a  second  creeping  wave-type  at 
f  -v.  572  MHz. 


For  exteided  dielectric  cylinders  in  a  lossy  medium,  the 
resonances  are  complicated  but  constructive  interference  of  re¬ 
flections  fron  the  front  and  rear  of  the  cylinder  yields 


cn 

4rVT 


♦ 
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which  for  an  ^  of  10  yields  (in  MHz) 

n  2  inch  3  inch  4  inch 

1  3164  2128  1580 

2  6328  4256  3160 


The  width  of  the  trenches  dug  to  bury  the  targets  was  1.33  feet  for 
the  dielectric  targets  and  2  feet  for  the  conducting  targets  and  the 
dummy  trench.  Reflections  from  the  side  walls  of  the  trench  would 
yield  resonances 


* 


clijtU  H2> 

2w/v 


where  w  is  the  width  of  the  trench, 
n  1.33  feet 


Again  for  an  er  of  10  (in  MHz) 
2  feet 


1  119 

2  357 


79 

237 


There  is  also  the  possibility  of  parallel  plate  waveguide- type 
resonances  associated  with  the  depth  of  the  trench.  This  effect  would 
produce  resonances  at  frequencies 
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where  d  is  the  trench  depth.  Resonances  at  multiples  of  33  MHz  for 
the  5  foot  depth  and  11  MHz  for  the  3  foot  depth  would  be  predicted. 
We  would  also  anticipate  seeing  at  least  the  first  resonance  of  the 
probe  structure  which  In  free  space  would  be  at  50  MHz.  An  exami¬ 
nation  of  the  various  spectra  shows  an  Invariant  resonance  at  roughly 
25  MHz,  and  this  resonance  Is  attributed  to  the  probe  structure  when 
it  is  on  the  ground.  One  possible  means  for  reducing  the  effect  of 
this  resonance  is  the  use  of  lossy  wire  so  that  the  signal  reflected 
from  the  end  of  the  probe  Is  attenuated  an  order  of  magnitude  before 
reaching  the  feed. 

Examining  the  spectra  in  terms  of  the  above  predicted  resonances, 
there  appears  to  be  a  correlation  between  the  trench  width  and  the 
resonances.  That  is,  the  spectra  corresponding  to  dielectric  targets 
all  have  a  resonance  near  119  MHz,  and  the  conducting  targets  near 
79  MHz.  We  would  not  anticipate  seeing  creeping  wave-type  resonances 
for  conducting  targets  with  the  power  and  target  depths  involved, 


54 


and  clearly  the  predicted  resonances  for  the  dielectric  targets  are 
beyond  our  capability  at  this  time.  There  Is  some  correlation 
between  the  predicted  and  measurod  dipole-type  resonances  for  the 
conduct  ng  targets,  e.g.,  the  1  foot  and  10  foot  (Figs.  23  and  26). 
»  VS,n^0nsistent<  The  sPectra  f°r  the  conducting  targets 
j;"!!™4*  but  the  differences  are  not  consistently 
related  to  the  difference  In  target  length. 


The  various  predicted  resonances  for  subsurface  cylindrical 
targets  can,  at  best,  be  considered  first  order  approximations  to 
more  complicated  mechanisms.  Other  processing  of  the  spectra  Is 

5wSibI?unn*2‘*  a.matcb?d  fHter  approach  using  predicted  resonances 
for  a  given  target  or  the  spectrum  of  the  dummy  trench.  Properly, 
however,  it  Is  felt  that  such  schemes  should  await  data  from  a  new 
set  of  measurements  where  higher  power  and  lower  frequencies  are 
available. 


tfll.  PULSE  SOUNDING  DATA  FOR  A  CYLINDRICAL  TARGET. 

INHOMOGENEOUS  STRATIFICATIONS  AND  A  CYLINDRICAL 
VOID  IN  A  LIMESTONE  MEDIUM 

earlJe.r*  the  original  research  plan  included  remote 
measurements  on  a  full  scale  basis  at  a  limestone  quarry*  located 

SSL5*?1*!!?  fr?I  the  Electr°Sc1ence  Laboratory.  Since  electrical 

tfiS  S!t?  w?f  not  available»  It  was  necessary  to  refurbish 

ThTJbIIhnl9f!^il!any  U5ed  for  terrain  c^ss  section  measurements. 8 

inlprlS^JSnh80!!  ty  cons1sting  of  a  2ls  ton  van  with  an  Independent 

nnfcl"?1  C3TbUS512n  gfnerator  Is  now  operational;  thus  electromagnetic 

LahorJnUndin9  data  at  a''ar1ety  of  sites  remote  from  the  ElectroScience 
Laboratory  are  now  possible. 

The  mobile  facility  has  been  used  to  obtain  initial  pulse  sounding 
data  for  particular  targets  in  a  limestone  medium.  Because  of  the 
!UI  Problems  discussed  earlier,  the  initial  measurements  were  made 
sing  the  identical  probe  structure  and  pulse  generator  utilized  for 
the  overburden  measurements.  A  photograph  of  the  mobile  rig  which 
houses  a  generator,  stabalines,  as  well  as  the  pulse  generator, 
sampling  scope  and  recording  equipment  is  shown  in  Fig.  29. 

mpHi.i!  fch°ct;9rapl?  sb°wing  the  typical  stratifications  of  the  limestone 
uhfrh^uac  sh?w"  30*  The  photograph  also  shows  a  large  tunnel 

which  was  originally  selected  as  one  of  the  geological  features  to 

J?te  tha*  tbe  mediuni  is  ^  no  means  solid  and  homogeneous, 
the  stratifications,  which  are  primarily  inhomogeneous,  are  typical  of 


♦Marble  Cliff  Quarries  Company,  2100  Tremont  Center, 
Columbus,  Ohio  43221 
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the  medium  in  the  locations  where  measurements  were  made.  In  Figs.  31a 
and  b,  the  response  of  the  medium  (no  target)  in  the  direct  reflection 
mode,  both  probes,  is  shown.  Note  that  the  response  waveform  is  lamer 
in  magnitude  by  rouqhly  a  factor  of  in  than  that  for  the  same  probe 
and  pulser  over  the  overburden.  The  mismatch  at  the  hnlun  ami  probe 
feed  point  is  larger,  and  a  large  reflection  from  the  end  of  the  probe 
is  clearly  evident.  It  is  obvious  from  these  waveforms  that  con¬ 
siderably  more  high  frequency  energy  is  arriving  at  the  end  of  the 
probe  than  was  the  case  for  an  overburden  medium.  Figure  31c  shows 
the  orthogonal  mode  response  in  the  same  location  (stratifications 
only).  Note  again  the  increase  in  the  magnitude  of  the  response. 

Figure  32  shows  the  same  sequence  of  response  waveforms  when  the 
probe  is  located  near  a  target  consisting  of  a  1  foot  diameter  hollow 
metal  pipe  whose  top  is  some  2  to  2.5  feet  below  the  limestone 
interface.  The  sketch  in  Fig.  32  shows  the  probe  and  pipe  con¬ 
figuration.  Some  evidence  of  the  target  is  apparent  if  these 
waveforms  are  compared  ti  the  no-target  set.  Finally,  in  Fig.  33, 
the  same  sequence  of  waveforms  when  the  probe  is  directly  over  the 
pipe  target  is  shown.  The  sketch  in  Fig.  33  shows  the  probe  and 
target  geometry.  Note  that  the  probe  is  positioned  for  maximum 
response  in  the  orthogonal  mode  and  one  would  not  anticipate  a  large 
target  response  in  the  direct  mode.  In  the  orthogonal  mode,  the 
presence  of  the  pipe  is  very  evident.  Figure  34  shows  the  amplitude 
spectra  of  the  cross  polarized  response  with  the  incident  pulse 
removed  for  a)  the  no-target  case  and  b)  over  the  pipe.  The  presence 
of  the  target  is  also  evident  from  a  comparison  of  the  S|  *ctra. 

Assuming  an  i  r  of  9  for  the  limestone  medium,  creeping  wave-type 
resonances  for  the  pipe  would  be  predicted  at  100  MHz,  218  MHz,  etc. 

The  100  MHz  resonance  can  be  seen  in  Fig.  34b.  The  same  sequence  of 
measurements  was  repeated  with  the  probe  located  above  the  tunnel 
shown  in  Fig.  30,  these  waveforms  are  shown  in  Fig.  35.  While  both 
the  direct  and  orthogonal  waveforms  indicate  responses  from  the 
homogeneous  and  inhomogeneous  stratifications,  there  is  no  evidence 
of  a  response  from  the  tunnel.  The  distances  are  such  that  a  response 
from  the  tunnel  would  occur  at  a  time  well  beyond  the  observed  duration 
of  the  waveform.  This  was  not  unexpected  considering  the  probe  and 
pulser  which  were  used.  The  limestone  quarry*  has  a  number  of  in¬ 
teresting  geological  features  suitable  for  tests  of  a  full  scale 
probe  and  a  pulse  generator  with  higher  power  and  lower  frequencies. 

An  extensive  series  of  measurements  at  these  locations  is  planned  for 
the  next  contract  period. 


*It  is  appropriate  to  acknowledge  the  excellent  cooperation  of 
Mr.  Max  Werner  and  his  staff  at  the  quarry  site.  Not  only  was 
free  access  to  the  site  provided,  but  Mr.  Werner  also  conducted 
tutorial-type  tours  of  the  quarry  for  our  education. 
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Fig.  32.  Response  waveforms  near  cylinder  target 
in  limestone  medium. 

(a)  direct  mode 

(b)  direct  mode 

(c)  orthogonal  mode. 
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Fig.  33.  Response  waveforms  over  cylinder  target  in  limestone  medium 

(a)  direct  mode 

(b)  direct  mode 

(c)  orthogonal  mode. 
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Fig.  35.  Response  waveforms  over  tunnel  in 
limestone  medium. 

(a)  direct  mode 

(b)  direct  mode 

(c)  orthogonal  mode. 
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Some  further  comments  concerning  data  processing  methods  for  this 
particular  electromagnetic  system  are  appropriate.  The  point  was 
made  earlier  that  a  proper  processing  must  await  an  analytical  capability 
for  a  specific  target  geometry,  i.e.,  a  calibration  capability.  With 
this  in  hand,  one  can  then  begin  to  secure  refined  estimates  of  the 
impulse  response  of  the  scatterer  from  measured  response  waveforms. 

Note  that  even  this  waveform  will  be  range  dependent,  not  only  because 
of  the  changing  characterise cs  of  the  medium  with  depth  but  because 
we  are  using  the  near-zone  fields  of  the  probe.  No  attempt  has  been 
made  here  to  optimize  the  processing  used.  Our  goal,  at  this 
stage,  was  simply  to  illustrate  the  system  capability.  Two  obvious 
things  could  be  done  to  further  demonstrate  the  detection  and  identi¬ 
fication  capability  of  the  system.  First,  a  bandpass  filtering  of 
these  same  data  based  on  the  predicted  resonances  given  earlier  is 
feasible.*  This  assumes  a  priori  knowledge  of  the  basic  target 
geometry.  Second,  a  much  more  refined  gating  of  the  target  and  no 
target  waveforms  or  the  difference  waveform  could  be  used.  With  this 
approach  a  measured  or  estimated  delay  time  for  the  medium  is  used 
with  the  gate  to  isolate  a  particular  portion  of  the  waveform.  This 
approach  has  been  successfully  used  by  Sullivan^  on  subsurface  targets 
at  shallow  depths.  An  a  priori  knowledge  of  the  target  geometry  is  not 
required.  It  was  not  felt  that  even  these  refinements  were  justified. 

A  few  conclusions  are  in  order  to  relate  the  system  discussed  in 
this  report  to  conventional  radar  systems.  The  major  difference  lies 
in  the  fact  that  present  studies  have  not  involved  targets  in  the  far 
field  of  the  antenna  and  that  the  limit  on  sensitivity  is  imposed  by 
ground  clutter,  not  the  signal  noise  level.  This  means  that  the 
received  power  level  does  not  vary  as  the  fourth  power  of  the  range  and 
also  that  the  ambient  noise  is  not  a  significant  factor.  These  factors 
will  indeed  become  important  as  the  range  is  increased.  The  fact  that 
clutter  is  the  factor  that  limits  sensitivity  also  implies  that  cor¬ 
relation  techniques  will  not  be  useful  in  improving  the  signal  to 
clutter  ratio  unless  the  characteristics  of  the  scatterer  and  the 
clutter  are  well  established.  This  is,  of  course,  our  goal  but  these 
data  at  this  time  are  not  available. 


*Successful  filter-type  processing  of  these  same  data  has  subsequently 
been  reported  JO 
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VIII.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  feasibility  of  an  electromagnetic  pulse  sounding  system  for 
the  detection  of  subsurface  targets  has  been  demonstrated,  in  the 
field,  from  measurements  of  conducting  and  dielectric  cylinders 
buried  in  the  overburden  and  from  measurements  of  a  conducting  cylinder 
buried  in  a  limestone  medium.  A  full  scale  version  of  a  probe  suitable 
for  use  with  the  system  has  been  designed  and  tested  with  low  power. 
These  tests  demonstrated  that  electromagnetic  pulse  energy  could  be 
effectively  coupled  into  a  lossy  medium  and  that  the  probe  was  very 
insensitive  to  obstructions  on  the  surface.  Because  of  delays  in  the 
delivery  of  special  high  power,  broad  band  baluns  needed  to  operate 
the  full  scale  probe  with  a  high  power  pulse  generator,  the  field 
measurements  were  made  with  a  scaled  down  version  of  the  probe  and 
a  low  power  pulse  generator.  The  pulser  output  power  and  spectrum 
were  inappropriate  for  interrogation  of  large  targets  at  deep  depths, 
nevertheless,  the  measurements  conclusively  demonstrate  the  feasibility 
of  the  pulse  sounding  approach. 

Analytical  studies  conducted  during  the  first  portion  of  the 
contract  have  developed  computational  tools  for  the  analysis  of 
various  probe  geometries  and  the  analysis  of  plane  wave  scattering 
from  planar  and  spherical  conductivity  contrasts.  The  comprehensive 
analysis  and  computer  programs  for  arbitrary  wire  antennas  in  a 
homogeneous  dissipative  medium  are  felt  to  be  a  major  advance  in 
the  state-of-the-art.  It  remains  to  add  a  half-space  correction  to 
these  programs. 

In  summary,  excellent  progress  has  been  made  toward  the  develop¬ 
ment  of  a  geological  tool  suitable  for  the  detection  and  delineation 
of  geological  or  man-made  anomalies  within  the  earth.  A  major 
demonstrated  feature  of  the  tool  is  its  ability  to  either  "see"  or 
be  "blind"  to  homogeneous  geological  features. 

It  is  recommended  that  during  the  next  contract  period  an  ex¬ 
tensive  series  of  measurements  be  made  in  a  local  limestone  quarry 
using  a  full  scale  probe  and  high  power  pulse  generators.  The 
necessary  equipment  for  these  tests,  including  a  mobile  rig  and 
special  high  power,  broad  band  baluns,  are  now  in  hand.  The  com¬ 
putational  programs  developed  from  our  analytical  studies  for  probe 
design  and  scattering  analysis  should  be  improved  and  extended.  In 
particular,  a  half-space  correction  for  the  wire  probe  analysis  and 
scattering  analysis  for  cylindrical  and  other  less  regular  targets 
are  needed.  Assuming  the  measurement  program  outlined  above  yields 
good  data,  more  advanced  processing  schemes  including  the  generation 
of  target-characteristic  response  waveforms  should  be  tested.  Finally, 
during  the  latter  stages  of  the  next  contract  period,  the  necessary 
plans  and  budget  estimates  to  permit  the  transport  of  our  equipment 
to  a  mine  or  tunnel  site  to  obtain  pulse  sounding  data  at  such  a 
location  should  be  made. 
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ABSTRACT 


A  concise  discussion  of  research  progress  on  Contract  H0210042 
for  the  period  22  February  1971  to  22  August  1971  is  presented. 

Analytical  results  on  the  scattering  by  various  planar  and  spherical  con- 
cuctivity  contrasts  and  on  design  data  for  an  electromagnetic  pulse 
sounding  probe  are  described  and  illustrated.  A  first  generation  version 
of  the  probe  is  given  and  initial  measured  data  demonstrating  certain 
features  of  the  probe  are  presented. 
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I. 


technical  report  summary 


detection  an^d^agnosTro^aeolooi^i SOUJdln9  te<*nique  for  the 
the  earth  is  being  investigated  9  ThP  U?nImade  anomalies  within 
technique  is  as  a 9h azard^ecti on  ton  T*  3*  applicatio"  of  this 
rapid  tunneling  operations  bStnJJ.  !  fiance  of  hard  rock 
visioned.  The  aoals  ni  ?h«  the^  p?tent1?1  uses  can  be  en- 
to  design  and  build  an  electroma^nPt?Unn9Kthlk-COntract  period  are 

couples  electromagnetic  eneJaTi2?o  JHCpPn  e  ®ffectively 
surface  targets  and  with  thiTn^h!  the  9round>  1s  insensitive  to 
from  cylinders  in  a  reck  2dlu S  reaSUre  the  scattered  fields 

ui  ti  e  s  ^th  at  tepresln  t 1  haza  rds  %o  r  ^ha  re/  rock  6f  ^  te  v 1 00  0f  discontin- 
The  effort  has  been  separated  into  S  r°Ck  tuf1neli ng  operations, 
the  other  experimental  While  it  ic  appr?acbes»  one  theoretical, 
«ork  Is  the  Jest  important  part  of  the  ^•'Z?d  Lhat  the  experiment 
gram  has  been  initiated  first  wh-tio  project,  the  theoretical  pro- 
interfaced  so  that X  eSeri  jSlli  SS^S-1"8"  obtained 
can  be  successfully  pursued^  ?futhls  ****rch  effort 

of  a  source  which  produces  a  nerindi^f"1  *hat_has  been  used  consists 
are  coupled  to  theqround  bv  J  brn^K  t!iain  of  vldeo  pulses-  These 
respond  to  the  ground  interface ^nd^hS^Mr3^'  The  probe  does  not 
cessfully  coupled  into  the  amnnd  nu  +kU  k6t^  ener9y  has  been  sue- 
contained  in  ?he  tratJtted^Sse  NoJ  does™??  band 

to  large  structures  above  the  around  I  !-U  appear  to  respond 
detected  when  driven  between  the  antenn?  fff  °?  wag°n  ha*  not  been 
currently  being  pursued  to  achievp  tho  c  ements.  Techniques  are 
medium.  The  probe  consists  nf  an  -h  !ame  9°al  for  a  rock  subsurfaa 
fed  dipoles  stret?heT^S  Sthe  suPr?;cre°off0tlh°SOn??,y  °riented  ba1“" 
launched  on  one  dipole  at  the  feed  ?  earth*  A  wave  ** 

ground.  As  it  propaqates  alnnn  +hnP2-  \tha*  penetrates  into  the 
trate  deeper  into  the  earth  ThphQhdlP°  uts  fringe  fields  pene- 
to  reduce  reflections  from  the  ante™?  s?f  t®  dl?ole  has  been  modifie 
term  nated  by  a  ground  ™re  a?  e  P  [U  fltUre  \tSelf*  lt  is  then 
on  the  wire  are  minimized  by  this  tlmS- 1  i-aJs  •  Sections 
all  of  the  energy  is  indeed  being  S  e  It  sl how;>  that  al™st 

transmitter  probe  and  the  orthoaona?  £™hp  *  Jhe  ground*  Both  the 
reflected  siqnals  Thp«:p  k  probe  can  be  used  to  receive 

acts  as  a  receiver.  In  our  presents vs temththSamplin9  Scope  wblcb 
the  sampled  voltage  and  stores  this  inf^f^®  50mputer  unitors 
and  can  plot  any  received  waveform  ThlTjtl?  f°,r  furtber  Processing 
non-homogeneous  targets  and  hence  is  idea?  £5  °9^  pr°be  "sees"  only 
cylindrical  obstruction.  ldeal  for  detect1on  of  any 
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Our  measurements  revealed  that  a  component  of  the  probe  system 
(balun)  which  matches  the  output  characteristics  of  the  pulse  generator 
to  the  input  characteristics  of  the  probe  was  not  capable  of  handling 

^cP0We;  f  generator  for  any  significant  pulse 

width.  This  made  the  actual  detection  of  subsurface  targets  with  this 
P^e  difficult.  This  is  not  a  basic  problem  because  other 
components  with  higher  power  handling  capability  are  available  and 
have  been  ordered.  A  different  feed  arrangement  which  eliminates  the 

aic  lhe  n*tch™Q  component,  albeit  with  a  somewhat  higher  mismatch, 
has  also  been  devised  In  the  runtime,  a  smaller  version  of  the  probe 
and  a  much  narrower  pulse  with  a  higher  frequency  operating  band  have 
been  used  to  demonstrate  the  feasibility  of  the  pulse  soundinq 
approach  with  measurements  of  a  clay  drain  tile  at  a  depth  of  3  feet 
This  probe  has  not  at  this  time  been  properly  matched.  However  by 
nreans  of  a  differencing  process  in  the  computer,  it  has  been  used  to 
ThfJ  ry  dlSt?nct  scattered  pulse  from  the  three  foot  deep  drain. 

TMs  differencmg  is  a  very  valuable  tool  which  can  be  retained  in 
the  final  system  by  a  simple  memory  and  a  simple  circuit  as  a  final 
toing  rechwism.  IhVec5uique  used  here  was  t0  "easure  the  apparent 
obtained  when  the  Probe  1S  Placed  over  a  section  ofground 
"h  C  ,therAre  no  Jar,9ets-  ^cse  reflections  from  the  probe  are 

tilP  Thla«?id  ^  th?se  wben.the  Probe  is  placed  over  the  drain 
le  The  reflected  pulse  obtained  in  this  manner  is  10  dB  above  the 

dUf  ,ThlsJnew  P^he  eliminates  the  need  of  grounding  the 

mHc  dlE0le  and  wou1d  be  adequate  for  penetration  of  the  hard 

rock  medium.  Preparations  for  measurements  at  a  nearby  site  which 

3nomaJies  characteristic  of  hard  rock  tunneling  hazards  have 

fnr  nnndtiaiT5-Se  prepa*a5u°ns  Primarily  involved  equipment  and  circuitry 

£  essentia  in^hl?^®  ??nputer’  The  versatility  of  the  computer 

to  fix  thl  ^  h  9  t-6  proper  Parameters  to  be  monitored  and 

°f0x  th®.  P^er  data  processing  steps,  i.e.,  design  of  the  system 

PvnonL  and  venfied  using  computer  software  without  the 

°  construction.  Once  the  proper  data  processing 

establ]sh®d>  tbese  should  be  replaced  by  devices  designed 
or  those  specific  tasks  in  the  final  system  for  hazard  detection. 

Preparations  are  also  being  made  at  the  nearby  quarry  for 
measurements  above  a  rocky  media  with  no  lossy  soil  layer.  This  has 
and*!?]  advanta9®J  and  disadvantages.  The  advantage  is  that  losses 

scatte^d  nu?JpeIf?f h  Wl1  -6  rdUCed  and  consequently  a  more  distinct 
thtrJho  1  *  1  be  ^c^ved  from  buried  objects.  The  disadvantage 

y?5e!"  can  not  be  as  easily  matched  and  above  ground 
reflection  could  be  more  apparent.  y 

-  The  feasibility  of  electromagnetic  pulse  sounding  as  a  tool 

earth  is1  also  °l  ge?logical  anomalies  within  the 

1S  *■  be  established.  A  pulse  sounding  approach  for  the 

i subsurface  targets  is  predicated  on  certain  advances 
in  identification  and  discrimination  of  radar  targets  and  on  advances 
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ih^ht!I!na  ^Sign  *ir  C0UP]in9  energy  into  a  lossy  medium;  both  of 

3J  tJ1s  1«borat°ny-  The  approach  also  utilizes 
c°Tter  technology  for  averaging,  differencing  and  other 
processing  of  periodic  signals.  In  a  field  version  system,  the 

nh^cot!h0functl?ns  wou)d  be  ^P^eeed  by  circuitry,  but  in  the  research 
p  ase  the  computer  again  offers  an  irreplacable  versatility  in 
processing  techniques.  y 

Simply  stated,  it  has  been  demonstrated  that  the  physical 

SriPmiltPi15  2Vn  objectJs!ze*  genenal  shape  and  composition)  are 
adequately  defined  by  the  interactions  of  the  object  with  electro- 

fS  S*2V?r  9  frec^ency  range  where  the  object  size  ranges 

KL5"11  tbat  comparable  to  the  interrogating  wavelength.  In 
terms  of  geological  anomalies  of  reasonable  size,  this  dictates  low 
frequencies  and  it  is  precisely  these  frequencies  which  can  penetrate 

fnr  th^intp  S1gnjficanb  dePths.  Thus  the  approach  is  ideally  suited 
for  the  interrogation  of  subsurface  targets, 

thp  c  ?Lth-S  stage  °T  the  contract,  analytical  tools  for  studying 
'  °J  electromagnetic  waves  by  planar  and  spherical 

approximations  of  geological  anomalies  and  for  the  analysis  of 

anni?r2mag?htlC  p^e  structures  have  been  developed  and  are  being 
applied.  The  need  for  these  calculations  is  to  establish  better  9 

fnr ^p/!Ia^ter-/0rJhe  sy$tem  and  to  provide  additional  data 
tor  target  identification  purposes.  The  latter  study,  essential lv 

°f  arb1trary  wire  structures  in  a  dissipative  medium/ 
is  felt  to  represent  a  significant  advance  in  the  state-of-the-art. 

forUintp2ener^0r  S'P'  214A^  With  pulse  characteristics  suitable 

pnrt  terr0^  10n  of  ?"°malles  of  the  order  of  10  to  100  feet  at 
depths  possibly  over  100  feet  has  been  purchased. 

nf  an  J iJ 5 IT*' ry ’  *be  an*lyt*cal  tools  for  the  design  and  analysis 

,amPantne'  cCt:°??9netlC  pul!e  soundin9  Probe  have  been  developed.  These 
same  tools  will  serve  to  define  the  basic  limitations  of  the  pulse 

sounding  technique.  A  first  version  of  the  probe  has  been  constructed 
a  power  handling  problem  revealed  and  solved,  and  the  feasibility ’ 
of  the  pulse  sounding  approach  demonstrated  via  measurements  of  a 
shallow  subsurface  target.  The  program  is  proceeding  in  agreement 

possibly  ars?inhl  Iffar-h  Eh3" .a?d  no. deviation  from  this  plan  except 
is  anticipated?*  th*  init1atlon  of  remote  site  measurements 


II.  PURPOSE 

.  ^-A  uni<?ue  electromagnetic  pulse  sounding  technique  for  the  detection 

hpindla9n°S^  9ec^°9'ical  and  man-made  anomalies  within  the  earth  is 

a  az9arndVdetPr9ta?ed\  ^  plication  of  this  technique  is  as 

azarci  detection  tool  in  rapid  excavation  operations  in  a  hard  rocP 

medium.  During  this  contract  period  the  objectives  of  the  pregram 
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are  to;  1)  design  and  demonstrate  a  probe  structure  capable  of 
effectively  coupling  electromagnetic  pulse  energy  into  the  ground  and 
at  the  same  time  be  insensitive  to  surface  scatterers,  2)  interrogate 
with  this  probe  certain  simple  subsurface  targets  such  as  voids  and 
discontinuities  representative  of  rapid  excavation  hazards,  and 
3)  establish  preliminary  bounds  on  the  workable  depths  of  electro¬ 
magnetic  pulse  sounding  via  theoretical  and  experimental  methods. 


III.  INTRODUCTION 

The  ultimate  goal  of  the  research  on  this  program  is  to  in¬ 
corporate  state-of-the-art  advances  developed  at  this  laboratory, 
both  in  antenna  systems  and  processing  techniques,  into  a  geological 
tool  capable  of  delineating  the  size  and  shape  of  geological  anomalies 
from  the  spectral  content  of  an  echo  pulse. 

If  the  gross  features,  i.e.,  overall  size  and  shape  of  an 
anomaly  are  to  be  deduced  via  electromagnetic  interrogation  then 
frequencies  over  an  8  or  10:1  bandwidth  with  a  fundamental  such  that 
the  linear  extent  of  the  anomaly  is  approximately  1/10  of  a  wave¬ 
length  in  the  medium  are  dictated. [1]  With  a  pulse  sounding  approach,* 
this  fixes  the  characteristics  (rise  time,  repetition  rate,  etc.)  of 
a  pulse  generator  suitable  for  a  given  range  of  anomaly  sizes.  It 
does  not  necessarily  follow  that  all  of  the  higher  spectral  content 
are  needed  to  obtain  a  viable  signature, [2]  but  the  lowest  fre¬ 
quencies  are  essential.  Of  course,  the  more  high  frequency  content 
available  the  sharper  the  reflected  pulses  and  consequently  an 
improved  range  gating.  This  is  an  advantage  of  the  hardrock  medium 
over  the  soil  media  currently  being  used  in  our  experimental  studies. 
Thus  the  approach  is  ideally  suited  for  subsurface  interrogation  where 
low  frequencies  are  needed  to  obtain  significant  penetration.  Clearly, 
the  attenuation  and  dispersion  of  the  ground  as  well  as  the  size  and 
scattering  characteristics  of  a  given  anomaly  bound  the  maximum 
working  depth  of  a  pulse  sounding  system.  Analytical  studies  to 
secure  estimates  of  these  limiting  effects  via  plane  wave  scattering 
computations  for  planar  and  spherical  conductivity  contrasts  are 
described  in  Section  IV. A  of  this  report. 

In  Section  IV. B,  a  comprehensive  analysis  of  wire  antennas  in 
a  homogeneous,  dissipative  medium  is  described.  Considerable  effort 
has  been  expended  on  this  study  since  realistic  estimates  of  pene¬ 
tration  depths,  dispersion  effects  and  surface  target  sensitivity  for 
practical  probe  geometries  are  a  vital  requirement  of  the  program. 

A  major  feature  of  the  computational  programs  developed  is  that  the 
speed  and  accuracy  are  sufficient  to  permit  pulse  type  excitation 
studies  via  Fourier  synthesis. 


*An  alternative  system  for  radar  applications  uses  discrete  harmoni¬ 
cally  related  cw  measurements  over  the  10:1  bandwidth. [3] 
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..  .  Cer*ain  exact,  closed  form  results  which  can  be  obtained  for 
the  transient  fields  of  an  infinite  line  source  in  the  pfesence  o? 

re?I  u  S"M"  Sec‘  °"  IV-C-  These  two-dimensional 

♦  IS  !  Erectly  applicable  to  the  iirmediate  problems  and 
have  not  therefore  been  extended  at  this  time.  P  ana 

thP  generat1on  of  the  pulse  sounding  probe  and 

These  in  tialft^tQU^mentSw0ri-thik  Pr°be  are  described  in  Section  V. 
nese  initial  tests  were  made  in  the  iirmediate  vicinity  of  the 

aboratory  but  the  necessary  circuitry  and  equipment  for  remote  site 
operation  are  now  nearly  completed.  ror  remote  si te 

our  0lJnJU!5r?h!f  the.present  status  of  this  research  effort  and 
and  Vi!  respectively™"”"9  C°"tract  period  are  gfven  1n  s«tions  VI 


IV.  ANALYTICAL  STUDIES 

A.  Scattering  Computations 

1 •  Planar  contrasts 

ho.  Th!  p1j”e  ”ave.scattering  characteristics  of  a  planar  contrast 
between  two  dissipative  media  are  most  simply  studied  in  terms  of  the 

transfLrSf0?he  ’Javeform  at  the  interface,  i.e.,  the  inverse  Laplace 
transform  of  the  frequency-dependent  Fresnel  reflection  coefficient 

For  normal  incidence  and  assuming  the  constitutive  parameters  of 

wa°  forermdlisas^owfnre?nUe^  the  «  tKTjSS  Ssponse 

waverorm  is  shown  in  Fig.  1.  The  response  to  any  other  incident 

wave'i r"i nci dent"^  “"S1*!0?  as  **  discussed*  ate/  The  p?ane 

of  an  impulse  singularity^?  height  medlUm  2'  Th®  reSP°"Se  consists 

(1)  A  = 


a  jump  discontinuity  of  magnitude 


(2) 


_ pz  !l_‘ 

‘of^TT +  LEr2  erl_ 


and  then  a  smooth  decay  to  zero  in  a  time  less  than  the  relaxation  time 
can^he  1®di urn  2.  Note  that  both  the  singularity  and  discontinuity 

to  ,0.0  with  the  °-’ 
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IMPULSE  RESPONSE  WAVEFORM  FOR  PLANE  WAVE  REFLECTION  AT 
THE  INTERFACE  OF  TWO  LOSSY  MEDIA.  NORMAL  INCIDENCE. 


*1  * 


A  -  WEIGHT  OF  IMPULSE  SINGULARITY,  A  -  frLlJ1** 


B  -  INITIAL  VALUE  OF  WAVEFORM  R  =  _ 

2 

ErlEr2 

ajl 

w  mi  i  lnu  vnuuL  u i  nnvuru i\pi «  o 

£o( 

Fri  +  Frt] 

r 

_er2  er2j 

T  IS  LESS  THAN  THE  RELAXATION  TIME^)  OF  THE  MEDIUM 


Fig.  1.  Sketch  of  Impulse  response  waveform 
of  planar  contrast. 
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Details  of  the  decay  rate  of  the  waveform  are  best  obtained  via 

t  uctTdVr^  :is:sip:ocedure,C4.]  wher*  the  "*«»"**  irin-8 

fr°  calculations  of  the  frequency-dependent  reflection 
coefficient  at  harmonically  related  frequencies  Examoles  of  curh 
waveforms  are  shown  in  Figs.  4  and  5.  To  si>?ify  the  synthesis 
procedure,  the  impulse  singularity  has  been  subtracted  from  the 

auleLtl0-,C°efficient  before  ^thesis.  The  synthesis  prScedur^ 

parameteTrin^nth0^!!^0^06  fre?u< ^-dependent  constitutive 

coefficient  2tbd^r^tdlafS1nCe  c?lculations  of  the  reflection 
coerricient  at  discrete  frequencies  are  made.  In  this  case  the  mn 

stitutive  ^rameters  used  to  calculate  the  quantities  in  Easfll 
limit  Sh°Uld  be  th°Se  the  redia  apProach  in  the  high  frequency 

*Pcn«nr1V|n  the  response  waveform  for  a  planar  contrast  the 

deDendenrI°r  an  1r)terrogating  plane  wave  with  an  arbitrary  time’ 
th?t  tho  o  ia!-  y  obtamed  by  convolution.[l]  Note  specifically 

pje^et0i„tch?^r^tved?:er 

aboIe  ^d?u™7  Si^o^drrco“c?eT!?etrSea21a,r'e?rth 

ation.  then  the  impulsive5  assonants  2  g25  one  °  *  tU""e,,n9  °per" 

diagnJJic7a?LCrLraS.KtS  *  d,stir":t,''e  signature  which  has 
sinaularitvfand  Jho  •  th®. S19n  and  magnitude  of  both  the  impulsive 
of  the  wlJefo™  ftJ?cP  J1SC0^tinuity  and  in  the  effective  duration 

1roulseWresponse  w^mls"  3  reaS°nab,e  app™(raatl'°"  for  the 
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As  with  the  planar  contrast,  the  assumption  of  a  plane  wave  incident 
on  the  target  places  restrictions  on  the  ambient  medium.  With  these 
programs,  formulation  of  the  transient  response  waveforms  via  Fourier 
synthesis  is  feasible.  At  present  the  computations  have  been  confined 
to  single  frequencies,  primarily  to  ensure  that  all  bugs  have  been 
removed  from  the  programs. 

Examples  of  the  computations  for  spherical  targets  are  shown 
in  Figs.  6,  7,  and  8.  In . Fig.  6  the  E-plane  (6a)  and  H-plane  (6b) 
patterns  showing  the  magnitude  of  the  scattered  field  squared  as  a 
function  of  bi static  angle  are  shown.  Zero  degrees  is  backscatter. 

The  sphere  in  Fig.  6  is  a  lossless  dielectric,  er  =  2.8,  with  a 
radius  of  5  cms.  The  ambient  medium  has  a  wavenumber  k|p  =  09.48  + 
jC.2973  so  that  the  index  of  refraction  for  the  sphere  is  1.056  - 
j0.003.  The  frequency  is  3.0  GH7.  Patterns  are  shown  in  Fig.  6 
for  observer  ranges  of  0.1,  0.3  and  0.5  meters.  In  Fig.  7,  similar 
results  for  a  lossy  sphere  (er  =  11.0,  a  =  0.385,  loss  tangent  = 

0.210)  with  a  radius  of  5  cms  are  shown.  The  ambient  medium  and 
the  frequency  are  the  same  as  those  in  Fig.  6.  Finally,  in  Fig.  8, 
results  are  shown  for  the  same  conditions  as  Fig.  7  but  a  lossy 
layer  1  cm  thick  (er  =  11.0,  o  =  0.0025,  loss  tangent  =  0.00136)  has 
been  added  to  the  sphere.  It  should  be  noted  that  for  economy 
reasons  the  angular  increment  for  these  computations  (Figs.  6,  7  and  8) 
was  quite  large  (30°),  thus  some  fine  detail  of  the  curves  may  not 
be  shown.  The  data  shown  were  primarily  run  to  check  out  the  com¬ 
puter  programs.  These  programs  are  efficient,  but  the  cost  of  such 
calculations  is  not  trivial.  Thus,  some  effort  to  estimate  parameters 
for  the  ambient  medium  and  spherical  targets  which  are  realistic 
models  of  tunneling  hazards  will  be  expended  before  extensive 
calculations  are  made. 

3.  Overburden  effect  estimates 

The  effects  of  the  air-earth  interface  and  propagation  in  the 
ground  on  the  plane  wave  interrogation  of  a  buried  target  has  been 
estimated  in  the  following  manner.  A  plane  incident  wave  at  an  off- 
normal  angle  from  the  earth  surface  is  assumed.  The  transmitted 
field  at  a  given  depth  is  calculated  under  the  assumption  of  frequency 
independent  constitutive  parameters  for  the  earth.  This  field  is 
then  modified  by  the  Mown  scattering  characteristics  (amplitude  and 
phase)  of  the  target  in  free  space  at  an  equivalent  (same  electrical 
size)  frequency.  The  propagation  path  to  the  surface  is  retraced  and 
the  field  then  transmitted  through  the  earth-air  interface.  It 
is  desired  to  compare  the  ramp  response  waveform  of  the  target  in 
free  space  and  in  the  eirth.  This  is  done  via  a  Fourier  synthesis 
procedure  with  calculations  as  described  above  at  discrete  harmonically 
related  frequencies.  For  the  ramp  response  waveform  (twice  integrated 
impulse  response)  only  a  very  few  frequencies  are  needed. [2]  It  is 
interesting  that  two  different  ramp  response  waveforms  for  the  target 
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.  Magnitude  squared  of  scattered  field  vs 
angle  for  a  lossy  dielectric  sphere  in 
medium,  (a)  E-plane,  (b)  H-plane. 
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can  be  obtained.  In  each  case  the  fundamental  frequency  is  calculated 
such  that  the  requirements  on  electrical  size  in  the  medium  described 
in  the  Introduction  are  met.  This  calculation  is  made  under  the 
assumption  that  the  earth  is  a  conductor,  i.e., 


(4)  x  medium  =  —  . 

/flfu  am 

o  m 

If  the  required  harmonic  wavelengths  are  calculated  from  Eq.  (4) 
then  one  response  waveform  results.  This  is  case  of  the  ramp 
waveform  incident  at  the  target.  If  one  simply  takes  multiples  of  the 
fundamental  then  a  different  waveform  is  obtained.  This  is  the  case 
of  a  ramp  waveform  on  the  transmitting  cable.  The  difference  is 
essentially  that  of  taking  the  harmonics  as  n^  or  n2^  (Eq.  (4)) 
where  <*>«  is  the  fundamental  angular  frequency.  In  practice,  either 
definition  could  be  used  provided  the  period  T  of  the  transmitted 
waveform  is  sufficiently  large  that  the  reflected  signal  from  the 
target  is  completely  damped  out  before  the  onset  of  the  next  pulse. 

This  implies  a  narrowly  spaced  spectrum  over  the  frequency  band  of 
interest  and  consequently,  one  could  reconstruct  via  the  computer 
either  of  the  incident  ramp  waveforms.  In  Figs.  9,  10  and  11  the 
free  space  and  underground  ramp  response  waveforms  for  a  conducting 
sphere,  conducting  cylinder  and  conducting  disk  are  shown.  The  time 
scale  in  these  figures  is  in  units  of  the  period  of  the  fundamental 
synthesis  frequency.  The  ground  parameters  are  taken  as  er  *  10.0, 
a  -  6.11  x  10“4  and  the  targets  are  at  a  depth  of  30  meters.  The 
interrogating  waveform  (periodic  with  ramp  discontinuities)  is 
incident  at  an  angle  30°  from  normal  with  the  electric  field  parallel 
to  the  ground.  The  targets  have  a  characteristic  dimension  of 
6  meters,  i.e.,  sphere,  disk  and  cylinder  have  a  radius  of  3  meters. 

The  cylinder  has  a  2:1  length  to  diameter  ratio.  The  axis  of  the 
cylinder  is  parallel  to  the  ground  interface  and  to  the  incident 
electric  field,  the  disk  is  perpendicular  to  the  ground  interface  with 
the  plane  of  the  disk  parallel  to  the  incident  electric  field.  10 
harmonic  frequencies  were  used  to  produce  the  waveforms  shown.  For 
each  target  it  is  seen  that  the  interface  and  ground  propagation  has 
distorted  the  free  space  waveform,  yet  the  distortion  is  not  such  as 
to  preclude  a  signature  assignment.  The  known  relationships!^]  between 
the  physical  properties  of  the  target  and  its  free  space  ramp  response 
waveform  have  been  distorted  and  clearly,  as  was  anticipated,  develop¬ 
ment  of  such  relationships  for  subsurface  targets  will  be  more 
difficult  than  for  the  free  space  case.  There  is,  however,  no  evidence 
to  indicate  that  the  effects  of  the  interface  and  overburden  will 
negate  a  time  domain  signature  approach. 
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Fig.  9. 
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(b)  underground  -  n^  harmonics, 

(c;  underground  -  n^  harmonics. 


1  Fig.  11.  Ramp  response  waveforms  of  a  conducting  disk: 

(a)  free  space, 

(b)  underground  -  nu0  harmonics, 

j  (c)  underground  -  n2u,0  harmonics. 
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B. 


Probe  Analysis  and  Design 


We  have  developed  a  comprehensive  analysis  and  computer  program 
for  wire  antennas  in  a  homogeneous  dissipative  medium.  The  program 
handles  all  types  of  wire  antennas  including  linear  dipoles,  V  dipoles, 
rectangular  loops,  circular  loops  and  arrays.  The  output  provides  a 
complete  analysis  of  the  system,  including  self  impedance,  mutual 
impedance,  current  distribution,  near- zone  fields,  and  far-field  patterns. 

The  finite  conductivity  of  the  wire  antenna  is  taken  into  account 
via  the  surface-impedance  formulation.  If  some  portions  of  the  antenna 
are  insulated  from  the  conducting  medium  by  a  thin  dielectric  sleeve, 
this  is  taken  into  account  via  the  equivalent  polarization  currents. 

A  piecewise-sinusoidal  expansion  is  employed  for  the  unknown 
current  distribution  on  the  antenna,  and  Galerkin's  method  is  used  to 
reduce  the  integral  equation  to  a  system  of  simultaneous  linear 
equations.  The  analysis  takes  place  in  the  frequency  domain.  The 
speed,  accuracy  and  generality  are  sufficient,  however,  to  permit  a 
Fourier  transform  to  the  time  domain. 

Figure  12  illustrates  typical  data  obtained  with  this  computer 
program.  This  figure  shows  the  near-zone  electric  field  distribution 
Ex  of  a  U-shaped  wire  antenna  at  100  KHz.  The  observation  point  is  on 
the  z  axis  w!iere  Ey  vanishes  and  Ez  is  much  smaller  than  Ex.  The  solid 
curve  shows  the  field  distribution  of  the  uninsulated  antenna,  and  the 
dashed  curve  shows  the  effect  of  insulating  the  horizontal  portion  of 
the  antenna.  The  input  power  was  one  watt  in  both  cases.  The 
polystyrene  insulation  has  an  inner  diameter  of  3/8  inch  and  an  outer 
diameter  of  3/4  inch. 

The  data  in  Fig.  12  indicate  that  the  uninsulated  antenna  may 
be  advantageous  for  detection  of  targets  at  close  range.  Of  greater 
importance  is  the  fact  that  the  computer  program  can  be  employed  to 
optimize  the  antenna  design  for  a  given  target  position,  frequency 
range,  and  soil  conditions. 

The  expansions  published  by  Banos  will  be  programed  to  account 
for  the  effects  of  the  air-earth  interface.  When  this  is  completed, 
the  computer  program  will  be  applied  to  the  analysis  and  design  of 
probes  for  the  electromagnetic  pulse-sounding  system. 


C.  Transient  Fields  of  a  Line  Source 

Thus  far  we  have  considered  plane  wave  incidence  for  the 
geometries  of  interest.  It  is  important  to  obtain  incident  fields  for 
different  sources  since  it  is  not  desirable  to  illuminate  the  target 
with  a  plane  wave.  While  the  technique  just  discussed  gives  the 
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fields  for  the  actual  sources  of  interest,  it  does  this  in  a  quanti¬ 
tative  way.  It  is  equally  important  to  evaluate  the  fields  of  various 
sources  in  functional  form  so  that  a  mote  thorough  understanding  of 
the  physical  mechanisms  involved  can  be  obtained.  The  line  source 
above  an  interface  is  most  important  since  simple  exact  closed  form 
expressions  can  be  obtained. 

„  Let  !  octangular  coordinate  frame  be  oriented  with  the  xy  plane 
on  the  surface  of  a  half-space  with  the  z  axis  perpendicular  to  the 
half-space  and  positive  into  the  half-space.  A  harmonic  line  source 
is  oriented  parallel  to  the  y  axis  at  a  height  h  above  the  half- 
space.  The  half-space  is  assumed  to  be  homogeneous  and  isotropic 

paran«te^  M.  MO.  O].  If  the  observation  point  is 
confined  to  the  z  axis  (x=0)  and  the  half-space  is  lossless  (oi=0) 

m6  Sh0Wn  •h-tJ  -han?e  0f  variable  and  Laplace  transformation 
of  the  well  known  seim-infimte  integral  solutions  for  the  fields[51 
yields  exact,  closed  form  expressions  for  the  real  time -dependent 
torment  fields  corresponding  to  a  step  current  (I  =  Inu(t))  excitation 
of  the  line  source.  The  y  component  of  the  electric  field  above  the 
half-space  (z<0)  is  given  by 


Jiwn  °f  the  electric  field  within  the  half-space  (z>0)  is 


For  the  same  excitation  and  a  lossless  half-space  an  exact  closed 
form  solution  can  also  be  obtained  when  both  the  observation  point 
and  line  source  lie  on  the  half-space  (z=h=0).  crvdtlon  P°int 
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where  i’(v,x)  is  the  incomplete  gamna  function 
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The  results  in  Eqs.  (8)  and  (9)  constitute  long  time  estimates  of  the 
transient  fields  for  a  step  current  excitation  of  the  line  source. 


If  the  half-space  is  homogeneous  with  frequency  independent 
parameters  then  the  response  waveforms  can  at  best  smoothly  decay 
from  the  short  time  estimates  in  Eqs.  (5)  and  (6)  to  the  long  time 
estimates  in  Eqs.  (8)  and  (9).  This  suggests  that  these  results 
might  be  combined  to  yield  a  simple  composite  waveform  which  could 
then  be  transformed  back  to  the  frequency  domain.  In  our  studies 
of  the  scattering  by  objects  in  free  space,  such  time  domain  com¬ 
binations  of  low  and  high  frequency  estimates  have  been  shown  to 
yield  surprisingly  good  approximations  in  mid- frequency  ranges.  It 
is  not  intended  to  further  pursue  this  topic  during  this  contract 
period,  but  such  an  approach  Would  appear  to  merit  reasonable 
conside  ation  in  a  longer  range  program. 


IV.  EXPERIMENTAL  STUDIES 

Near  the  beginning  of  the  contract  period,  preparations  for 
experimental  measurements  at  a  remote  site  (a  nearby  quarry)  which 
offered  hazard-type  anomalies  were  initiated.  A  pulse  generator 
(H.P.  214A)  appropriate  fo'1  such  measurements  was  ordered.  This 
pulser  has  a  peak  pulse  amplitude  of  50  volts  into  50  ohms  or  100 
volts  in  1200  ohms,  a  15  ns  rise  time,  pulse  widths  variable  from 
45  ns  to  the  ^s  range  and  repetition  rates  from  10  kc  to  1000  kc. 

The  major  preparations  involved  replacing  a  direct  interfacing  of  the 
antenna  system  with  an  instrumentation  computer,  which  had  been  used 
on  earlier  studies,  with  a  recording  arrangement  such  that  a  non-real 
time  use  of  the  computer  was  feasible.  It  should  be  understood  that 
in  an  ultimate  system  the  computer  functions  will  be  replaced  by 
circuitry.  At  this  stage,  however,  the  computer  offers  an  irreplacable 
flexibility  in  processing  changes  and  experimentation.  The  above 
mentioned  preparations  are  now  nearly  completed  -  certain  components 
have  not  yet  been  delivered.  The  system  for  remote  recording  of  data 
resamples  the  oscilloscope  output  waveforms  at  a  very  low  speed  with 
computer  control  of  the  sampling  point  relative  to  a  trigger  synched 
with  the  scope  sweep.  As  shown  in  Fig.  13  the  trigger  and  sampling 
point  offset  voltage  are  fed  into  a  phase  locked  loop.  The  voltage 
controlled  oscillator  (VCO)  output  drives  a  one  shot  which  provides 
the  proper  pulse  shape  for  triggering  the  sample  and  hold  amplifier. 

The  offset  volti.-,;  varies  the  phase  of  the  VCO  output  relative  to 
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Fig.  13.  Data  conversion  circuit  for  tape  recorded 
remote  site  measurements. 


the  trigger,  thus  effectively  moving  the  sampling  point  across  the 
signal  waveform.  The  sampled  signal  is  fed  to  the  computer  analog/ 
digital  input  for  conversion  to  digital  form. 

A  breadboard  version  of  the  device  is  currently  being  built. 
?h!mh1e  1ine!rV?te9rat?d  circuit  Performs  the  entire  function  of 
and  hold6"  °Cked  °°P*  A  mimature  hybn'd  TOdule  Performs  the  sampl 

in  thoThfL?5Cill°LC°P!  waveforms  wil1  be  recorded  on  magnetic  tape 
n  the  field  and  the  tapes  later  played  back  into  this  device  in  the 

IfJrihi’h  Sh"k“  5h1s  d?vice  does  "dt  hare  t0  9°  the  remie 

«  L  ^adb°ard  "i™0?  can  be  used  for  actual  data  conversion 
as  soon  as  it  is  operational. 

fh„  PU1se  generator  was  delivered  during  the  fifth  month  of 

and  S1"Ct  that  date  t6StS  0f  various  antenna  configuration: 

;nberr°gatl°"  of  certain  subsurface  targets  in  the 
vicinity  of  the  laboratory  have  been  underway  using  the  direct 
interfacing  to  the  instrumentation  computer. 

A.  First  Generation  Probe 

centerTJpHbHi,ioif1iC-r0ma9n!Jic  pr°be  being  studied  is  an  insulated, 
oilh!  ?d  dl?  1ying  on  the  ground  W1’th  the  ends  of  the  dipole 
either  directly  grounded  by  rods  driven  into  the  earth  or  capacitivlv 
oupled  using  plates  on  the  ground.  The  arms  of  the  dipole  may  be 
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single  wires  with  ground  rods  or  two  wires  in  a  U  or  V  arrangement  to 
two  separate  ground  rods.  Alternatively,  the  dipole  arm  may  be  a 
v-shaped  conductor  eliminating  the  need  for  a  ground  rod.  The 
system  presently  envisioned  consists  of  two  such  dipoles  arrayed 
orthogonally.  Two  modes  of  operation  are  possible;  either  direct 
reflection  using  one  or  the  other  of  the  dipoles  as  transmitter  and 
receiver  or  orthogonally  using  one  dipole  as  a  transmitter  and  the 
other  as  a  receiver.  Both  modes  are  important  since  the  direct 
mode  is  sensitive  to  symmetrical  planar  stratifications  whereas  the 
orthogonal  mode  is  blind  to  such  targets. 


A  sketch  of  one  version  of  the  probe  and  a  block  diagram  of  the 
instrumentation  is  shown  in  Fig.  14.  In  Fig.  15,  the  time  domain 
reflection  of  one  of  the  dipoles  shown  in  Fig.  14  is  shown  in  the 
direct  reflection  mode.  The  incident  pulse  for  the  results  shown  in 
Fig.  15  was  very  long,  essentially  a  step  with  a  15  ns  rise  time. 

Shown  in  Fig.  15  are  the  differences  when  the  arms  of  the  dipole  are 
a  single  wire  or  two  wires  in  a  U  or  V  arrangement.  Also  shown  are 
the  reflected  waveforms  when  the  feed  cable  is  terminated  in  an 
open  circuit,  short  circuit  and  a  matched  load.  By  adjusting  the 
ground  rod  depth,  the  latter  portion  of  the  waveform  can  be  reasonably 
matched  to  the  matched  load  waveform  indicating  that  at  low  frequencies 
at  least  the  antenna  is  well  matched.  The  effects  of  tuninq  by 
varying  the  ground  rod  depths  is  shown  in  Fig.  16,  these  results  also 
indicate  that  energy  is  being  coupled  into  the  ground. 

...  The  w?vef°rns  ™  Fi?s*  15  and  16  correspond  to  an  incident  pulse 

WJ^0a?nr2?l!I'ate]t!  V°lt  peak  amplitude*  Attempts  to  operate  the 
probe  in  either  the  direct  or  orthogonal  mode  using  the  peak  output 

of  the  pulse  generator  revealed  that  the  balun  (transformer)  used 
to  go  from  50  ohms  unbalanced  to  200  ohms  balanced  was  saturating  at 
a  leye1  we!!  beiow  the  peak  output  of  the  pulser.  This  is  not  a  basic 
difficulty  since  baluns  with  a  higher  power  handling  capability  are 
available  and  have  been  ordered.  In  the  interim  before  their  delivery 
however,  the  pulser  cannot  be  used  at  output  levels  appropriate  for 
interrogation  of  subsurface  targets.  Tests  did  establish  however  that 
the  probe  was  insensitive  to  surface  structures  (an  automobile  driven 
into  various  quadrants  of  the  probe  had  no  effect  for  a  low  level 
incident  pulse)  and  that  a  synthetic  target  in  the  form  of  two 
conducting  plates  connected  electrically  could  be  detected 


The  temporary  balun  difficulty  was  circimvented  by  qoinq  to 
essentially  a  smaller  version  of  the  probe  (the  dipole  arms  became 
solid  conductors  in  a  V  shape)  and  the  use  of  a  different  pulse 
generator.  The  pulser  is  not  appropriate  for  deep  penetration  but 
can  be  used  to  interrogate  relatively  shallow  targets.  Usinq  this 
system,  a  clay  drain  tile  buried  at  a  depth  of  3  feet  has  been 
detected.  The  equipment  consisted  of  an  pulse  generator,  a  "four- 
leaf  clover  probe  lying  on  the  ground  surface,  and  a  computer- 
controlled  sampling  oscilloscope.  The  pulse  was  approximately 
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46»?  BLOCK  DIAGRAM  —  ORTHOGONAL  MODE 


Fig*  14 •  (j)  Electromagnetic  probe. 

(b)  Block  diagram  of  pulse  sounding  system. 
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Fig.  15.  Direct  reflection  mode  of  one  probe 
various  terminations. 


Fig.  16.  Direct  reflection  mode  of  one  probe  -  effect  of  ground  rod  depth. 


2  nsec  long  with  a  peak  amplitude  of  10  volts,  and  a  pulse  repetition 
rate  of  approximately  1  MHz.  A  sketch  of  the  antenna  is  shown  in 
Fig.  17. 

The  antenna  was  located  in  several  positions  in  the  vicinity 
of  a  storm  drain-tile  line  running  behind  the  laboratory  building 
at  about  a  3  foot  depth.  Several  locations  over  the  drain  tile  and 
several  locations  of  assumed  "target-free"  ground  were  measured,  and 
the  differences  between  the  "target"  and  "no-target"  waveforms  wc re 
computed  and  plotted  as  shown  in  Fig.  18.  Difference  curves  for 
target-free  ground  are  shown  in  Fig.  19.  These  can  be  considered  to 
represent  the  clutter  level.  These  curves  demonstrate  that  the  pipe 
is  visible  to  the  system,  and  that  several  features  of  the  return 
waveforms  remain  constant  from  one  location  over  the  pipe  to  the 
next.  Based  on  the  variation  of  these  return  waveforms  it  is  esti¬ 
mated  that  the  signal -to-cl utter  ratio  of  the  pipe  signature  is  about 
10  dB. 


Several  measured  difference  waveforms  of  the  catch  basin  at  the 
end  of  the  drain  tile  are  shown  in  Fig.  20.  Comparison  of  these 
waveforms  with  those  of  Fig.  18  indicates  that  the  nature  of  the 
return  signal  is  sensitive  to  the  shape  of  the  underground  target. 
Thus,  identification  of  underground  object  shape  from  these  signature 
waveforms  seems  feasible. 


B.  Propagation  Test 

While  it  is  not  presently  possible  to  use  the  probe  in  Fig.  14 
with  high  power  because  of  the  balun  saturation,  the  dipoles  can  be 
used  as  receivers.  For  the  transmitter,  the  balun  requirement  can 
be  eliminated  by  going  to  an  unbalanced  antenna.  To  this  end  a 
monopole  antenna  consisting  of  a  2  foot  rod  on  a  circular  (1-1/2  foot 
diameter)  ground  plane  was  constructed.  In  operation,  the  ground 
plane  rests  on  the  ground  with  the  rod  driven  into  the  ground.  This 
antenna  can  be  fed  directly  by  coaxial  cable.  Results  of  transmission 
tests  rsing  the  monopole  as  a  transmitter  and  the  probe  of  Fig.  14 
as  a  receiver  are  shown  in  Fig.  21.  The  monopole  was  near  broadside 
to  one  dipole  and  near  endfire  to  the  other  and  located  at  a  distance 
of  160  feet  from  the  vertex  of  the  probe.  The  transmitted  pulse  had  a 
pulse  amplitude  of  50  volts,  a  pulse  width  of  50  nanoseconds,  and  a 
repetitive  rate  of  10  KHz.  It  was  established  that  the  received  wave¬ 
form  shown  was  not  the  effects  of  a  surface  wave  by  moving  a  large 
conducting  screen  in  the  vicinity  of  the  monopole  without  disturbing 
the  received  waveform.  It  is  not  yet  understood  why  the  received 
waveforms  are  identical  and  further  measurements  are  in  progress.  It 
has  not  been  established  as  yet  what  penetration  depths  are  being 
achieved;  with  the  monopole  as  a  transmitter  it  is  possible  that  energy 
is  being  confined  close  to  the  surface  of  the  ground.  Measurements  are 
presently  in  progress  to  establish  the  penetration  achieved  as  well 
a*  Che  attenuation  vs  frequency  for  the  local  ground. 


Catch  basin  target  difference  waveforms. 


Fi9-  21.  Received  pulse  waveform. 
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'I.  CONCLUSIONS 


The  plane  wave  scattering  by  planar  conductivity  contrasts  has 
been  explored  via  the  impulse  response  waveform  for  the  contrast. 

A  Fourier  synthesis  technique  is  utilized  to  determine  the  decay 
rate  of  the  waveform  and  will  be  used  to  explore  the  effects  of 
frequency-dependent  constitutive  parameters.  The  response  waveform 
has  a  distinctive  signature  and  diagnostic  features  for  both  the 
dielectric  constant  and  conductivity  contrasts. 

Computer  programs  for  computing  the  plane  wave  scattering  by 
arbitrary  spherical  targets  in  a  dissipative  mediun  have  been  written 
and  verified,  including  the  case  of  layered  spherical  targets.  Calcu¬ 
lations  for  realistic  spherical -type  tunneling  hazards  and  for  deducing 
characteristic  response  waveforms  for  such  hazards  are  in  progress. 

The  attenuation  and  dispersion  effects  of  a  lossy  overburden 
have  been  estimated  by  assuming  plane  wave  propagation  and  utilizing 
the  free  space  scattering  characteristics  of  spherical,  cylindrical 
and  disk  targets.  These  results  indicate  that  interface  and  over¬ 
burden  effects  do  not  preclude  unique  signatures  for  the  targets. 

Exact,  closed  form  expressions  for  the  transient  fields  of  a 
line  source  in  the  presence  of  a  dielectric  half-space  have  been 
obtained  for  particular  observer  locations.  Quasi -static  estimates 
of  the  transient  fields  for  a  lossy  half-space  have  also  been  obtained. 
From  these  results,  simplified  estimates  of  the  fields  of  a  line 
source  in  the  presence  of  an  arbitrary  half-space  appear  to  be 
feasible. 

A  comprehensive  analysis  and  computer  program  for  all  types  of 
wire  probes  in  a  homogeneous  dissipative  medium  have  been  developed. 
With  these  tools,  design  data  for  pulse  sounding  probes  to  interrogate 
at  various  depths  are  being  prepared. 

Initial  measurements  have  been  made  on  two  versions  of  the 
electromagnetic  pulse  sounding  probe.  The  results  indicate  that  the 
probe  can  be  well  matched,  does  couple  energy  into  the  ground  and  is 
insensitive  to  surface  targets.  The  presence  of  a  clay  drain  tile 
at  a  depth  of  3  feet  has  been  detected  using  a  small  version  of  the 
probe  and  a  pulse  generator  inappropriate  (at  least  where  a  lossy 
overburden  exists)  for  deep  penetration.  The  limited  power  handling 
capability  of  transformers  used  to  match  the  unbalanced  coaxial  feed 
to  the  balanced  probe  has  temporarily  precluded  use  of  the  full  scale 
probe  with  the  high  power  pulse  generator.  Larger  transformers  have 
been  ordered  and  other  schemes  for  eliminating  the  need  for  trans¬ 
formers  are  being  studied. 


Transmission  tests  demonstrating  that  pulsed  electtomagnetic 
energy  can  be  propagated  through  the  ground  over  considerable  distances 
(160  feet)  have  been  made.  The  energy  may  be  confined  to  shallow  depths, 
and  this  possibility  is  being  explored. 


VII.  FUTURE  PLANS 

It  is  intended  that  major  emphasis  during  the  remaining  contract 
period  will  be  concentrated  on  experimental  measurements  of  subsurface 
targets,  first  at  this  laboratory  and  then  at  a  remote  test  site 
where  anomalies  characteristic  of  tunneling  hazards  are  available. 

This  is  in  accordance  with  the  original  research  plan. 

Analytical  studies  will  be  primarily  devoted  to  the  wire  antenna 
analysis  from  which  realistic  dispersion,  attenuation  and  penetration 
calculations  as  well  as  probe  design  data  are  possible.  The  scattering 
from  spherical  contrasts  will  be  calculated  for  hazard-type  parameters 
and  response  waveforms  obtained  via  Fourier  synthesis.  Introduction 
of  frequency-dependent  constitutive  parameters  will  be  made  for  both 
the  antenna  analysis  and  spherical  scatterer  program. 
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